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Heat semigroup and singular PDEs 


I. BAILLEUlEand F. BERNICOT0 
(with an Appendix by F. Bernicot D. Frey) 


Abstract. We provide in this work a semigroup approach to the study of singular PDEs, in the 
line of the paracontrolled approach developed recently by Gubinelli, Imkeller and Perkowski. 
Starting from a heat semigroup, we develop a functional calculus and introduce a paraproduct 
based on the semigroup, for which commutator estimates and Schauder estimates are proved, 
together with their paracontrolled extensions. This machinery allows us to investigate singular 
PDEs in potentially unbounded Riemannian manifolds under mild geometric conditions. As an 
illustration, we study the generalized parabolic Anderson model equation and prove, under mild 
geometric conditions, its well-posed character in Holders spaces, in small time on a potentially 
unbounded 2-dimensional Riemannian manifold, for an equation driven by a weighted noise, 
and for all times for the linear parabolic Anderson model equation in 2-dimensional unbounded 
manifolds. This machinery can be extended to an even more singular setting and deal with 
Sobolev scales of spaces rather than Holder spaces. 
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Resolution of PAM in such a 2-dimensional setting 


Introduction 

1.1 State of the art Following the recent breakthrough of Hairer [35] and Gubinelli, 

Imkeller, Perkowski |32|, there has been recently a tremendous ac¬ 
tivity in the study of parabolic singular partial differential equations (PDFs), such as 
the KPZ equation 

{dt -dl)u = + ?, 

the stochastic quantization equation 

{dt - A)u = -u^ + ^, 

or the generalized Parabolic Anderson Model equation 

(dt - A)u = F(u)^ 

in all of which ^ stands for a space or space-time white noise. Each of these equations 
involves, under the form of a product, a term which does not make sense a priori, given 
the expected regularity of the solution in terms of the regularity of the noise Hairer’s 
theory of regularity structures is built on the insights of earlier works [MlIlQlETI on 
(1 -|- l)-dimensional space-time problems where he used the framework of rough paths 
theory, under the form of Gubinelli’s controlled paths, to make sense of previously ill- 
posed singular PDEs and give a meaningful solution theory. Rough paths theory was 
used in this approach as a framework for studying the properties in the 1-dimensional 
space variable of potential solutions. However, the very notion of a rough path is 
intimately linked with the 1-dimensional time axis that parametrizes paths. 

To by-pass this barrier, both the theory of regularity structures and the paracon- 
trolled approach developed in [32| take as a departure point the fact that, like in rough 
paths theory, to make sense of the equation, one needs to enrich the noise ^ into a finite 
collection of objects/distributions, and that one should try and describe the potential 
solution of a singular PDE in terms of that enriched noise. The latter depends on the 
equation under study and plays in the theory of regularity structures the role plaid by 
polynomials in the usual world to give local descriptions of functions under the form 
of Taylor expansions at every space-time point. The description of a solution in the 
paracontrolled approach is of a different nature and rests on a global comparison with 
the solution to a linear equation, (dt — A'ju = ^, in the above examples, via the use of 
Bony’s paraproduct. In both approaches, the use of an ansatz for the solution space 
allows for fixed point arguments to give a robust solution theory where the solution 
becomes a continuous function of all the parameters of the problem. 

So far, both theories have only been formulated and tested on some singular PDEs 
on the torus, to the exception of the works |38l [39] on the parabolic Anderson model 
equation in and , and our follow up work [7] . We introduce in the present work a 
functional analytic setting in which we are able to extend the paracontrolled approach 
of [32] to investigate singular PDEs of the form 

(dt + L)u = F{u,^) 

for a second order differential operator L, and a nonlinear term F(tt,^), on potentially 
unbounded, Riemannian or even sub-Riemannian, manifolds or graphs. (The change of 
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sign — to + in the operator is irrelevant.) This is a priori far from obvious as the main 
analytic tools used in the paracontrolled approach in the torus involve technics from 
Fourier analysis that do not make sense on manifolds or graphs. We develop to that 
end a functional calculus adapted to the heat semigroup associated with the operator 
[dt + T) , which we use to define a paraproduct enjoying the same regularity properties 
as its Euclidean analogue. Such paraproducts adapted to a semigroup, as well as a par- 
alinearization theory, have already been studied in recent works [iniiii]. However, the 
irregular character of the noises ^ involved in the above motivating equations requires 
us to improve the definition of such paraproducts so as to build a framework where to 
consider regularity with a negative exponent; such an extension will be provided here. 
Building on these tools, one can set up, as in [32], a framework where to investigate 
the well-posed character of a whole class of parabolic singular PDEs. It is especially 
nice that all the objects in our framework are defined uniquely in terms of semigroups, 
unlike the notions of Holder spaces used in the theory of regularity structures that 
involve a metric structure unrelated to the equation under study. As a by-product, 
we are able to handle some general classes of operators L whose treatment seem to be 
beyond the present-day scope of the theory of regularity structures, as illustrated in 
some examples given in section 12.11 

It is unclear presently how one can adapt the different notions and tools of the theory 
of regularity structures to extend them to a (Lipschitz) manifold or graph setting, or 
to other second order operators (other than the Laplace operator), or how to work 
with Sobolev spaces instead of Holders spaces. Apart from the very definition of a 
regularity structure on a manifold, the existence of the reconstruction operator in this 
setting seems in particular challenging, as its proof in involves some deep results on 
wavelets that were not proved so far to hold true on generic manifolds, not even on all 
open sets of M'^. Their extension to a non-smooth setting also seems higly non-trivial. 
So it comes as a good news that one can use some reasonably elaborate theory of 
semigroups to implement the alternative machinery of paracontrolled calculus in that 
setting; as described below, it also allows us to have much flexibility on the operator L 
and also on the geometry of the ambiant space. Roughly speaking, we could say that the 
point of view of the theory of regularity structures relies on the metric and differential 
properties of the underlying space, while the present extension of the paracontrolled 
calculus corresponds to a functional point of view adapted to the operator L involved 
in the parabolic singular PDEs. We link here these two sides of the medal by requiring 
from the heat semigroup (e“*'^)f>o to have a kernel together with its gradient, that 
satisfies pointwise Gaussian bounds; this describes in some sense the link between the 
functional calculus and the ambiant space, with its metric and its differential geometry. 

We explain in Appendix B how this approach can be used in the context of Sobolev 
spaces rather than Holder spaces. The former setting is slightly more difficult to handle, 
from a technical point of view, since Sobolev spaces involve simultaneously all the 
frequencies, whereas for Holder spaces we can work at a fixed frequency scale. We do 
not know how such extension could be implemented within the setting of the regularity 
structure. 

The first part of this work is devoted to a precise study of the so-called paracontrolled 
calculus in a very abstract setting, given by a doubling ambiant space, equipped with 
a self-adjoint operator —L generating a semigroup with Gaussian bounds for its kernel 
and its gradient. A suitable definition of paraproducts is given, and the main rules of 
calculus for paracontrolled distributions are described. This general theory is all we 
need to study a number of parabolic singular PDEs on manifolds. 
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1.2 Paracontrolled calculus The mechanics of paracontrolled calculus [32] is elementary 

and easy to use; we describe it here as the present work only 
extends its scope while keeping its structure untouched - unlike our work [7|. We use 
a somewhat informal style in this section and take the stochastic PDE given by the 
parabolic Anderson model equation (PAM) 

(1.1) {dt + A)u = u^ 

as an illustration. The symbol ^ stands here for a 2-dimensional spatial white noise, 
of Holder regularity —1“. At the beginning of this story is the fact that one expects a 
solution u to that equation to be (—l“-|-2)-parabolic Holder regular, as a consequence of 
the regularizing properties of the heat semigroup, while this regularity is not sufficient 
to make sense of the product u^, as the sum of their parabolic Holder regularity is not 
positive. The point, however, is that u is not expected to be any kind of l“-Holder 
function, rather it is expected to behave, at small space-time scale, like the solution X 
of the elementary well-posed equation 

{dt + A)X = ^, 

with null initial condition. The paracontrolled approach to solving the 2-dimensional 
(PAM) equation consists in the following three step process. 

(1) Set yourself a .^-dependent ansatz for the solution space, made up of func¬ 
tions that behave like X at small space-time scale, and equipped with a 
Banach space structure. 

(2) Show that the product is well-defined for any element u of the ansatz 
space. 

(3) Solve the equation via a fixed point argument. 

The subtle point here is that the data of the distribution ^ itself is not sufficient to 
give sense to the product u^, and that we really need that ^ be random to build on the 
probability space where it is defined another distribution together with which one can 
make sense of all the products u^, for any u in the ansatz space. Once this enrichment 
of ^ has been done by purely probabilistic means, the above three step process is run 
in a deterministic setting. 

From a technical point of view, a l“-Holder function u will be said to behave like X 
at small scale if it is of the form u ~ n^(A), with v bounded, up to some term more 
regular than X] write u = n^(Ai) + u^, with a remainder of Holder regularity strictly 
greater than 1. The bilinear operator H. (•), which will be a generalized paraproduct, 
that appears here has good continuity properties on large classes of distributions and 
satisfies the identity 

ab = na(6) -I- nfe(a) -|- n(a, b) 

for any bounded functions a, 6, for a continuous operator n(-,-) on x that 
happens to extend continuously to pairs of Holder regular distributions for which the 
sum of their regularity is positive. In the torus, the Littlewood-Paley decomposition of 
a and h as an infinite sum of smooth functions whose Fourier transforms have supports 
in dyadic annuli can be used to define H.(•) and n(-, •), by writing 

ab = '^aibj= ^ aibj + ^ bjai+ ^ aibj. 
i,j i<j-l i<i-l 

This definition justifies that we call n(-,-) the diagonal operator. The following for¬ 
mal analogy with the rules of stochastic calculus will enlighten the core technical tool 
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of paracontrolled calculus described in a second. Recall that if M and N are two 
continuous martingales one has 

d{MN) = MdN + NdM + d{M, N). 

The above space of functions u = n^(X) + can be turned into a Banach space. 
Once this ansatz for the solution space has been chosen, remark that the product 
can formally be written as 

= n„(^) + n^(u) + n(u,o 
= njo + + n(n,(x), e) + n(u«,0- 

Since has Holder regularity strictly bigger than 1 and ^ is — l“-regular, the sum 
of their regularity indices is positive, and the term n(u^,^) is perfectly well-defined. 
This lives us with n(n„(X),,^) as the only undefined term. The following fact is the 
workhorse of paracontrolled calculus. The trilinear map 

C{a, b, c) := n(na(6), c) — a n(6, c) 

happens to depend continuously on a, b and c provided they are Holder distributions, with 
the sum of their Holder exponents positive. Note the paralell between the continuity 
of this ’commutator’ and the rule for stochastic differentials, for which, given another 
continuous martingale P, we have 

d(^j^ MdN ,P'^ = Md{N,P). 

The formal product can thus be written as a sum of well-defined terms plus the 
formal product un(X, ^), with a diagonal term n(X, still undefined on a purely 
analytic basis. This is where probability comes into play. If one regularizes ^ into 
with defined accordingly, one can prove that there exists a function/constant 
such that the renormalized quantity := converges in probability 

to some limit distribution of Holder regularity 0“ = 1“ this is enough 

to make sense of the product on an analytical basis; but replacing 

by ^(X^e) in the decomposition of amounts to looking at the product 

u — C^) . The enhancement ^ := of f is called a rough, or enhanced, 

distribution, and one can use it to define the product from the above formulas. 
At that point, it does not come as a surprise that one can then set (PAM) equation as 
a fixed point problem in the ansatz space, and that the unique solution to the problem 
(as it happens to be) is the limit of the solutions to the elementary problems 

{dt -h A)u^ = u^ (^^ - C^). 

The (PAM) equation (11.11) is said to have been renormalized. 

More complicated problems are treated along these lines of reasoning in ISailHllMl 
[65], to cite but a few. All use the above paraproduct machinery in the setting of the 
torus, where it is defined via Paley-Littlewood decomposition, such as described in this 
section. We introduce in this work a far more flexible paraproduct, defined intrinsically 
in terms of the semigroup associated with the operator L that plays the role of A in the 
equation, in a very general geometrical setting. This offers the possibility to investigate 
stochastic PDEs in a manifold setting, which is our primary motivation. We gain much 
flexibility along the way, in terms of operators that can be used in place of A, and even 
in Euclidean domains, the scope of the present work seems to be beyond the present 
day knowledge provided by the theory of regularity structures. 
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1.3 A generalized parabolic Anderson model As an illustration of our machinery, we 

study the stochastic PDE given by the 
generalized parabolic Anderson model equation (gPAM) 

dtu + Lu = F{u) n(0) = uq, 

on some possibly unbounded 2-dimensional Riemannian manifold M satisfying some 
mild geometric conditions. One can take as operator L the Laplace-Beltrami operator 
or some sub-elliptic diffusion operator; see section 12.11 for examples. The nonlinearity 
F is , and ^ stands here for a weighted Gaussian noise with weight in n L°° - see 
the definition in section 15.31 The deterministic counterpart of the (gPAM) equation 
can be set once we are given a rough distribution ( = (Cj we write informally 
for the product operation in the ansatz space - a slightly different and more precise 
notation will be adopted later on. The following results involve some parabolic Holder 
spaces , with negative exponents 7 , that are defined in section [ 23 ] in terms only of 
the semigroup Pt generated by — L. We refer to section [5] for a full statement and proof 
of these results. The geometrical assumptions stated here are introduced and explained 
in the next sections. The next two statements hold true provided the heat kernel of the 
semigroup generated by — L, together with its gradient, satisfy some Gaussian bounds 
given in the begining of section (Til conditions on L ensuring that these bounds hold are 
also given there. The letter C" stands for a spatial (that is time-independent) Holder 
space. 

► Theorem 1.1. Assume that the measured metric manifold (M, d, /i) is equipped with a 

volume doubling measure, and that the heat semigroup generated by —L has a kernel that 
satisfies some Gaussian bounds (UE), together with its gradient (Lip). Let (ii). 
an initial data uq € a nonlinearity F £ C^, and a positive time horizon T. Let 
C = be a rough distribution, with ( £ and 

(a) Local well-posedness for (gPAM). For a small enough time horizon T, the 
generalized PAM equation 

(1.2) dtu + Lu = F{u) ■ (, u{0) = Uq 

has a unique solution. 

(b) Global well-posedness for (PAM). Under the assumption that the rough 
distribution takes values in some space of weighted distributions, the PAM equation 

dtu + Lu = u ■ C, u(0) = uo 

has a unique global in time solution in some function space. 

The implementation of this result in the case where ^ ^ is a random Gaussian 

spatial noise takes the following form, for a precise version of which we refer to theorem 
15.51 it holds in the same geometrical setting as the above result. The property of the 
measure // put forward in the statement is called Ahlfors regularity; write Bj.{x) for 
a metric ball of center x and radius r. 

► Theorem 1.2. Suppose, in addition to the Gaussian bounds (UE) and (Lip) satisfied by 

the heat kernel and its derivative, that the reference measure /i on the manifold M is 
doubling and satisfies the uniform lower bound fj,(^Br{x)) > cir^, for all x £ M, for some 
positive constant ci and the homogeneous dimension v. Let ^ stand for a time-independent 
weighted noise in space, and set := and A^(f) = Pt-s(^^) ds. 
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(a) There exists a time-independent function C"^(-) := E n(L on M such 


that the pair —converges in probability to a random rough distribution 


(b) If stands for the solution of the renormalized equation 


dtu^ + Lu^ = F(u^) — C^¥'{u^) F(ti^), = uq 

it converges in probability to the solution u of equation (|1.2I) driven by 


Note that one cannot expect the renormalizing function to be constant unless 
the manifold M is homogeneous and the operator L commutes with the group action 
- which holds in the torus when working with the Laplacian. Note also that we do 
not assume M to be bounded. Working with a weighted noise rather than with white 
noise allows us to by-pass the somewhat heavy use of weighted Holder spaces, such as 
done in [38l|39] and [7]; the latter work deals, among other things, with paracontrolled 
calculus in weighted Holder spaces. 

We have organized our work as follows. Section [2] presents the functional setting 
in which our theory is set. The main geometrical assumptions on the geometric back¬ 
ground are given in section 12.11 where examples are given; these assumptions involve 
the properties of the heat kernel of the semigroup generated by L. A family 

of operators is introduced in section 12.21 which will play in the sequel the role played 
by Fourier projectors in the classical Littlewood-Paley theory. We introduce in section 
E3]a scale of Holder spaces, defined uniquely in terms of the semigroup (e 
A paraproduct is introduced in section 13.11 and is shown in section 13.21 to enjoy Fhe 
same continuity properties as its Euclidean analogue. A crucial commutator estimate 
between paraproduct and resonant terms is proved in section 13.31 together with some 
paralinearization and composition estimates in section [331 Following [32], we then in¬ 
troduce in section ITT] what plays the role in our setting of paracontrolled distributions, 
and prove some fundamental Schauder estimates in section 13.51 Sections 2 to 4 give us 
all the material needed to investigate a number of singular PDEs on manifolds from 
the point of view of paracontrolled distributions. Section [5] is dedicated to the proof of 
theorems 11.11 and 11.21 

We end this work by Appendix B, jointly written with Dorothee Erey, in which we 
explain how we can weaken our assumption of Lipschitz regularity of the heat kernel 
(Lip), which we make in the main body of this work, in terms of more geometrical 
properties. We also show that one can prove results in Sobolev spaces similar to those 
proved in Holder spaces in the main body of that work. 

We collect here a number of notations that will be used throughout that work. 

• Eor a ball B of radius r and a real A > 0, denote by \B the ball concentric with 
B and with radius Ar. We shall use u < v to say that there exists a constant C 
(independent of the important parameters) such that u < Cv and u ~ u to say 
that u < V and v < u. We also adopt the non-conventional notation 7 ^ for the 
classical gamma function, dehned for a > 0 by the formula 





the capital letter F will be used to denote the carre du champ operator of some 
other operator. 

• For p G [l,oo] and every / G L^, the L^-norm, with respect to the measure p, 
is denoted by ||/||p. For p,q G [l,oo], we write ||r||p^q for the operator norm 
of an operator T from to 

• For an integer k > 0, we write for the set of functions continuously differ¬ 
entiable /c-times / : M —)■ R, equipped with the norm 


ct 


oo + sup 




Functional calculus adapted to the heat semigroup 

As announced in the introduction, this section is dedicated to describing the func¬ 
tional framework where we shall set our study. Section 12.11 sets the geometrical frame¬ 
work needed for what we want to do, in terms of a semigroup. We introduce in section 
12.21 some operators that will play the role of ’localizers’ in frequency space. These oper¬ 
ators are used in section [2Al to define a scale of Holder spaces which will be instrumental 
in the sequel. 

2.1 Heat semigroup on a doubling space Let denote by {M,d) be a locally compact 

separable metrisable space, equipped with a 
Radon measure p, strictly positive on any non-empty open set. Given a ball B{x,r) 
of center x and radius r, the notation V{x,r) will stand in the sequel for p(^B{x,r)). 
To make things concrete, the space (M, d) will mainly be for us smooth Riemannian 
manifold or a (possibly infinite) metric graph. We shall assume that the metric measure 
space (M, d, p) satisfies the following volume doubling property 

(VD) R(x,2r)<R(x,r), 

for all X G M and positive r, which can be stated equivalently under the form 
(2.1) ^(a:,r)<(^^) V{x,s), 

for some positive scaling factor v, for all x G M, and all 0 < s < r; it implies the 
inequality 

for any two points x, y in M and 0 < s < r. (Another easy consequence of the volume 
doubling property is that balls with a non-empty intersection and comparable radii 
have comparable measures.) 

Let also be given a non-negative self-adjoint operator L on L‘^{M,p) with dense 
domain 'D 2 {L) C L‘^{M,p). Denote by £ its associated quadratic form, defined by the 
formula 

£{f,g) ■= [ fLgdp, 

Jm 

on a domain B which contains T> 2 {L). We shall assume that the Dirichlet form £ is 
strongly local and regular; we refer the reader to the books |26l IM] of Fukushima & 
CO. and Gyrya-Saloff-Coste for precise definitions and background on Dirichlet forms. 
(The reader unfamiliar with this setting may think of the Laplace operator in a com¬ 
pact Riemannian manifold.) These two properties will be obviously satisfied in the 
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examples we shall work with. It follows from these conditions that the operator L gen¬ 
erates a strongly continuous semigroup of contractions on L^(M,/i) which is 

conservative, in the sense that = 1, for all t > 0; see e.g. Subsection 2.2.7 in the 

book [31]. We shall also assume that the semigroup has a kernel, given for 

all positive times t by a non-negative measurable real-valued function pt on M x M, 
such that 

(e“*^/)(x) = / pt{x,y)f{y)dfi{y), 

Jm 

for /x-almost all x in M, and every / G 'D2{L). The kernel pt is called the heat kernel 
associated with L. We assume that it satisfies for all 0 < t < 1 and /x-almost all x, y, 
the following typical upper estimates 

Pt{x,y) < 

Jv{x,Vi)V{y,Vt) 


Under the volume doubling condition (|VDD . the previous estimate self-improves into a 
Gaussian upper estimate (UE) for the heat kernel and its time derivatives 

d{x,yf' 


(UE) 


dtPt{x,y) 


< 

rs_/ 


^Jv{x,^/t)V{y,^/t) 


exp —c 


that holds for a fixed positive constant c, for all integers a, all times 0 < t < 1, and p- 
almost every x,y £ M; see for instance the article |28l Theorem 1.1] for the Riemannian 
case, and the work [201 Section 4.2] for a metric measure space setting. We also assume 
that the heat kernel satisfies the following Lipschitz condition 


(Lip) 


Pt{x,y) -ptiz,y) 


< 


/ d{x,z) \ 

V Vt J 


1 

^Jv{x,Vt)V{y,Vt) 



Let insist here that inequalities (jUEp and (tnH) are assumed to hold only for 0 < 
t < 1, rather than for all positive times. It follows classically from the Gaussian 
estimates (jUEp and the volume doubling property that the heat semigroup 
is uniformly bounded on LP{M,p) for every p G [l,oo], and strongly continuous for 
p G [l,oo). Last, note that is, under these conditions, bounded analytic on 

L^{M,p), for every I < p < -|-oo, which means in particular that the time-derivatives 
are bounded on LP{M,p) uniformly in 0 < t < 1, for every integer 

n > 0; see [^ . 


A comment is in order here, about our two assumptions (lUEp and (Lip). In the theory 
of regularity structures or Euclidean theory of paracontrolled calculus, regularity at any 
order may be considered because of the implicit use of the very nice differential geometry 
of Euclidean space, or the torus. In our current and far more general framework, since 
we only have a pointwise assumption on the heat kernel and its gradient, it is natural to 
expect that one cannot quantify the regularity of some objects to an order greater than 
I. That is why in the different statements proved in the next sections some extra mild 
conditions on the regularity exponents will appear, as compared with their Euclidean 
analogue. Since we aim to work within the present optimal / minimal setting, these 
new limitations cannot be removed without additional assumptions, and we shall be 
restricted to study regularity properties at order at most I, including negative orders; 
this is not restrictive as far as applications are concerned in the present work. 


Here are four representative classes of examples of doubling metric measure spaces 
and Dirichlet forms satisfying the above conditions. This list of examples emphasizes 
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that we have much flexibility in terms of the operator L as well as in terms of the 
underlying space {M,d,n). 

(a) Markov chains. Let X be a countable set equipped with a Markov chain, 
specified by a symmetric Markov kernel k : X x X ^ ]R_|_, and let m be a 
non-negative function on X, used to define a measure m on X, with density m 
with respect to the counting measure fi. Denote by (•, •)m the scalar product on 
Consider also for integers n > 1 the iterated kernel k'^ defined recursively 
by k^{x,y) := f k^~^(x, z)k(z,y) fi{dz). Denoting by K the symmetric Markov 
operator with kernel k - with respect to y, the formula 


x^y^X 

fx — 

rrix 


kxyi^fx fy){9x 9y) 


= E 

xex 


9x 


^ ^ kxy9y j 
y&X 


mm 


= </> 

associated with the non-negative self-adjoint operator 

{L9){x) = ^{9x-Y 1 ^-y9y) =^{9x- (^<?) J, 

* y&X ^ 

defines a (strongly local) regular Dirichlet form and allows us to generate the 
continuous heat semigroup (The above sum in x is implicitly re¬ 

stricted to those X for which > 0, so there is no loss of generality in assuming 
that m > 0. ) The map k induces a distance d on X by letting be equal to 

{n > 1 ] 3 zq, ..., with zq = x, Zn = y and k{^Zi, ^i+i) > 0, for i = 0..n — l}, 

for y ^ X. Following Grigor’yan’s result [29], one can give growth conditions on 
the m-volume of d-balls that ensure the conservative character of the semigroup 
generated by L in £^(m). Then it is classical that getting Gaussian upper 
estimates for the semigroup is very closely related to getting discrete¬ 

time versions of Gaussian estimates for the iterated Markov chains 
and similarly for the Lipschitz regularity of their kernels. Usually, given such a 
discrete framework, one prefers to work with the discrete-time Markov chains 
rather than the continuous heat semigroup. To obtain upper Gaussian estimates 
and a Lipschitz regularity for the iterated Markov chains on a graphs is the 
topic of a huge literature to which we refer the reader; see for instance works by 
Hebisch and Saloff-Goste m for discrete groups and by Ischiwata [H] for an 
extension to nilpotent covering graphs and more recently [l5| for a perturbation 
of these previous results. For example, the regular graphs and (Z/XZ)'^ 
have heat semigroups satisfying the Gaussian estimates (lUED and the Lipschitz 


property ( Lip ). Needless to say, for a (large) finite graph {X,E), with edge set 
E, and b^y = 1 if {x,y) G E, and rux = Yly^x^xy, the previous results hold 
with the graph distance in the role of d. 


(b) Second order differential operators on Riemannian manifolds. Let (M, d, y) 

be a doubling, possibly non-compact,” complete Riemannian manifold with 
Ricci curvature bounded from below. Then the heat semigroup gen¬ 

erated by the Riemannian Laplace operator satisfies both the upper Gaussian 
estimates (|UEp and the Lipschitz regularity (Lip) for small time 0 < t < 1, and 
for every time f > 0 if the Ricci curvature is nonnegative; see [61] and [38| for 
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references. Particular examples are smooth compact Riemannian manifolds, or 
unbounded Riemannian manifolds with pinched negative Ricci curvature, such 
as hyperbolic spaces. 

Even on the Euclidean space we may consider a second order divergence 
form operator L = —div(AV) given by a map A taking values in real symmet¬ 
ric matrices and satisfying the usual ellipticity/accretivity condition. Then if 
A is Holder continuous, it is known that —L generates a self-adjoint semigroup 
satisfying the properties (IUE |1 and ; see [1]. Similarly, consider an open 
bounded subset H C (with Lipschitz boundary for example to ensure the 
Ahlfors regularity) and consider the self-adjoint Laplace operator L associated 
with Dirichlet or Neumann boundary conditions. There is an extensive litera¬ 
ture to describe assumptions on H such that (jUEIi and ( |LipD are satisfied. The 
present setting may well be beyond the present scope of regularity structures, 
for which the Green function of the operator needs to satisfy some regularity 
assumptions that were not proved to hold true under a sole Holder continuity 
assumption for A, and whose formulation on a manifold is a real problem out¬ 
side the realm of Lie groups or homogeneous spaces. On the other hand, the 
theory developed here works well in that relatively minimal setting. 


The estimates (jUEl) and (Lip) also hold when working on a convex or C^- 
regular bounded subset of the Euclidean space, with L given by Laplace operator 
with Neumann boundary conditions; see m- 


(c) Sub-elliptic left invariant diffusions on groups. Let G be a unimodular con¬ 
nected Lie group, endowed with its left-right Haar measure /r. Consider a family 
X := {Xi, of left-invariant vector helds on G satisfying Hormander con¬ 

dition. They dehne a class of admissible paths 7 ,, characterized by the existence, 
for each of them, of measurable functions oi, ..., 0 ^ such that one has 


7'(f) = ^ai(f)Xi(^(t)). 


2=1 


The length of such a curve is defined as 

u / i- 




and the (Carnot-Caratheodory) distance d{x, y) between any two points x, y of 
G is defined as the infimum of the lengths of all admissible curves joining x to 
y. We then consider the sublaplacian A defined by 




2 = 1 

Then the operator A generates a heat semigroup satisfying both the upper 
Gaussian estimates (jUEp and the Lipschitz regularity ( |LipP for small time 
t G (0,1]; see for instance Chapter 8 in the book [62]. If the group is nilpo- 
tent then it is also globally doubling [33] and so the heat semigroup satisfies 
the Gaussian upper bound (lUEp and enjoys the Lipschitz property ( Lip ) for 
every f > 0; see jOUllSS]. Particular examples of such groups, are stratified Lie 
groups, and so Heisenberg groups. Eor such Heisenberg-type Lie groups, a kind 
of Eourier transform may be defined involving their irreducible unitary repre¬ 
sentations, which can be used to define an analog of the Euclidean paraproducts 
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/ paradifferential calculus, such as done is m- We shall see, as a by-product of 
the present work, that the structure of heat semigroup is sufficient to construct 
similar tools with greater scope. 


(d) The general case given by a subelliptic operator is more difficult. Let (M, d, fi) 
be a complete and smooth connected manifold endowed with a self-adjoint 
smooth locally subelliptic diffusion operator L satisfying LI = 0. Then Bau- 
doin and Garofalo introduced in [S] a property, called “a generalized curvature- 
dimension inequality”, which has to be thought of as a lower bound on a sub- 
Riemannian generalization of the Ricci tensor. Under such a condition, the 
heat kernel generated by L satisfies (lUEp as well as ( |Lip ); see [M] • We refer the 
reader to [9] for some examples of such sub-elliptic settings and the fact that 
the heat kernel also satisfies in that case some Gaussian lower bound. 


Throughout that work, a point o £ M will be fixed, which we shall use to define a 
class of test functions, together with its ’dual’ class of distributions. 


► Definition. We define a Frechet space of test functions setting 


So:={f& f|p2(L'^);Vai,a2GN, (l + d(o, 


n>0 


< oo 


with 


ll/ll := sup lA 

ai,a2SN 


(l + d(o,-))“^L“V 


A distribution is a continuous linear functional on So', we write S'o for the set of all 
distributions. 


(We point out that the arbitrary choice of point o ^ M is only relevant in the case 
of a unbounded ambiant space M; even in that case, the set So does not depend on 
o, for o ranging inside a bounded subset of M.) Every bounded function defines for 
instance an element of 5'. Examples of test functions are provided by the pt(x, •), for 
every fixed x G M and 0 < t < 1. Indeed for integers 01,02, the upper bound (lUEp 
with the analyticity of the semigroup yield that •) satisfies the same upper 

Gaussian estimates than the heat kernel itself and so we deduce that 


(1 -F d{o, y)) •)) (y) 


<-—^(l + d(o,y)re-^^ 
V{x,Vt) 

f— 0.2 

< - 

V{x,Vt) 


/ d(x,y)^ 


for some positive constants c and cL Note that the heat semigroup acts not only on 
functions, but also on distributions, by setting 


((e^"<^),/) := (^,(e‘V)) 


for 4> G So and f G So- We refer the reader to m and [59] for more details on the 
extension of the semigroup to distributions. 

For a linear operator T acting from So to 5', it will be useful below, to denote by 
Kt its Schwarz kernel, characterized by the identity 


{T{f),g) = j KTix,y)f{y)g{x) y{dy)y{dx), 


giving an integral representation for every f,gG So- 
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2.2 Time derivatives and Carre du champ of the Let us introduce here a family of op- 
semigroup erators that will play the role in our 

setting of the Fourier multipliers used 
in the classical Littlewood-Paley theory, that localize a function in frequency space. 
These will be the building blocks used to define a convenient paraproduct for our 
needs, such as done below in section 13.11 


► Definition. Given a fixed positive integer o, set 
( 2 . 2 ) := 
and 

(2.3) := (j)a{tL), where := — [ —, x > 0, 

Ta Jx 

for every t > 0. 

So we have for instance = e~^^, and = tLe~^^. The two families of 
operators and are defined so as to have the relation 

(2.4) tdtP^^ = tL4>'a{tL) = - 7 -' qS“\ 

so = (—and = pa{tL)e~^^, for some polynomial pa of degree 
0 — 1, with Pa(0) = 1. The analyticity of the semigroup provides a direct control of the 
operators and 


Proposition 2.1. For any integer a > 0, the operators and Q'f'^ have kernels satisfying 
the Gaussian estimate (lUEp . and the Lipschitz regularity property ( |Lip ); as a consequence, 
they are bounded in every spaces for p G [l,oo], uniformly with respect to t G (0,1]. 


ha) 


Following the above interpretation of the operators and the following 

Calderon reproducing formula provides a decomposition of a function / in U’{M,p) 
into a low frequency part and a high frequency part very similar to the Littlewood-Paley 
decomposition of a distribution in terms of frequencies; see e.g. [8]. 


► Proposition 2.2 (Calderon reproducing formula). Given p G (1,-Foo) and / G L^{M,p), 
we have 


hmP/“V = / mLP{M,p) 

for every positive integer a, and so 


(2.5) 




Proof — One knows from theorem 3.1 in [23], that the operator L has a bounded 

functional calculus in LP{M,p) under the volume doubling condition on {M,d,p), 
and the assumption that the heat kernel satishes the upper estimate (lUEl) . Since 
this implies in particular sectoriality of L in LP{M, p), Theorem 3.8 in [21] yields the 
decomposition of LP{M, p) into nullspace and range of L. Using this decomposition, 
the Convergence Lemma implies for every / G U’{M,p) 

/ = hrnP/“V = - [" dtP}'^^fdt + Pi'^\f) 

Jo 

= 7a^ £Qi"Jj + Pi'"\f), 
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where the limit is taken in LP{M, and where we have used identity (12.4^ : see e.g. 
[31 Theorem D] or [361 Lemma 9.13]. 

I> 


We shall also make an extensive use in the sequel of the square-root of L, given by 
its carre du champ operator T, defined for all {f,g) G T>2(L) x 'D 2 {L) as a bilinear 
operator satisfying the identity 

£{f,g)-= [ fL{g)dn= [ gL{f)dfi= f T{f,g)dn. 

Jm Jm Jm 

It is also given by the explicit formula 

T{f,g) = -]^(LUg)-fL{g)-gL{f))- 

we shall write T>2(r) C for its domain, which contains T) 2 {L). 

As a shorthand, we write r(/) for r(/,/)2 in the sequel, which can be thought as 
the length of the intrinsic gradient of /. It follows from the conservative property of 
L and its non-negative character, that the bilinear map T is positive and satisfies the 
identity 

l|r(/)||'2=/ r{f,f)dfi= [ fL{f)dg = £{f,f). 

Jm Jm 

From its positive property, a generalized Cauchy-Schwarz inequality yields that for 
every f,g^ T>2(L) then 

(2.6) |r(/,5)P < r{f,f)T{g,g) = T{f)T{g). 

According to the Beurling-Deny-Le Jan formula, the carre du champ satisfies a Leib¬ 
niz rule 

(2.7) r{fg,h) = fr{g,h)+gT{f,h), 
for all f,g,h G 'D 2 {T), and a chain rule 

(2.8) L{F{f)) = F'{f) L{f) + F''{f) r(/, /). 

for every function F G CKM) and every / G 'D 2 {L); the function F{f) is automatically 
in V 2 {L) - see e.g. [26l Section 3.2] and [58l Appendix] for these points. 

The following pointwise and LP-estimate for the intrinsic gradient of the semigroup 
will be used several times in a crucial way; its proof is given in Appendix^ It says that 
the carre du champ of the semigroup satisfies also some Gaussian pointwise estimates, 
as given by the following claim. 


Proposition 2.3. The following inequality holds 

(2.9) 

(\/tr)(e"*^/)(xo) 


< 


^ Jv{xo,Vi)V{y,Vi) 


exp —c 


d{xo,yf 

t 


\f{y)\dy{y), 


for every t > 0, every function f G L'^, and almost every xq G M. Consequently, we have 


sup 

t>o 


(VJr)(e-‘V) 


p^p 
-tL 


< oo. 


for every p G [1, oo]. We may replace the semigroup e in the above equations by any of 
the operators Pf°'\ for any a > 0. 











15 


2.3 


Holder and Besov spaces through the heat 
semigroup 

space if 

II/IIa- := ll/lloo + sup 

0<d(x,y)<l 


Let us recall as a start that given a pa¬ 
rameter cj G (0,1], a bounded function 
/ G L°° is said to belong to the Holder 

- u -^- < OO- 

d{x,yf 


Recall on the other hand the definition of the inhomogeneous Besov spaces associated 
to a semigroup; they were precisely studied in several works, such as m or m , to 
name but a few. We shall make an extensive use of these spaces. 


Definition 2.4. Fix a positive integer a, an exponent p,q G (l,oo), and a G 
distribution / G 5', is said to belong to the Besov space ^ if 


B‘^ ■= e 

p,q " 


(i 






9 dt\ « 

— < oo. 

p t ) 


This definition of the space does not depend on the integer a > 1, provided a is big enough. 


We refer the reader to m for details about such spaces and a proof of the fact that 
they do not depend on the parameter a used to dehne them, provided a is sufficiently 
large with respect to a. The limiting case p = q = oo leads to the following definition. 


Definition. Let a positive integer a be given. For a G (—oo,2), a distribution / G 5' is 
said to belong to the space if 






-F sup 
oo 0<t<l 




t 2 < oo. 


This definition of the space does not depend on the integer a > 1. 


We give in Appendix [A] a simple and self-contained proof that the space does not 
depend on a, and that any two norms || • Wc’, defined with two different values of a, are 
equivalent. The following proposition justifies that we call the spaces Holder space, 
for all CJ < 2, possibly non-positive. 


► Proposition 2.5. For a G (0,1), the spaces and are the same and the two corre¬ 
sponding norms are equivalent. 


We give here a complete proof of this proposition as it provides an elementary illustra¬ 
tion of how the properties of the operators are used to make actual computations. 
This kind of reasoning and computations will be used repeatedly in the sequel, when 

working with our paraproduct. Recall that the operators have kernels K {a) sat- 

Qt 

isfying Gaussian pointwise estimates, by proposition 12.11 


Proof — We divide the proof in two steps, by showing successively that A'^ is continu¬ 
ously injected in , and that, conversely, is continuously injected in . 

Step 1. ^ . Note first that since the Holder space A”" is made up of bounded 

functions, it is included in S'^. Fix an integer a > 1; then for every t G (0,1), we 
have 

(Qi“V)(a:) = - f{x))Yx) = J K^(a){x,z){f{z) - f{x)) p{dz). 

For the points z G M, with d{x, z) < y/i < 1, we have 

\f{z) - f{x)\ < d{x,zyW/W a^ < t^ll/IU- 
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SO that 


/ K (a) {x,z){f{z) - f{x)) fi{dz) <t2||/||A-/ K (a) {x,z) fx{dz) 

Jd{x,z)<\/t * J ^ 


since has a kernel satisfying Gaussian pointwise bounds. The same bounds 
show that 

K (a) {x,z) {f{z) - f{x))n{dz)\ 




f \/t<d{x,z)<l 


< IIJIIA' 


I y/t<d{x,z)<l 


J V 

<i^ll/llA- 

Similarly, we have 

/ “ f{x))d'{dz) 

Jl<d{x,z) * 


y/t<d{x^z)<l ^(3^7 V^) 


K^{a){x,z) d{x,zY ^{dz) 

1 




< 


< 


' l<d(x,z) 


K^(a){x,z) fi{dz) 


e ‘ J A- 


so it comes that the inequality 

(qS“V )(3 


<t2 


< t 2 




A'^ 


holds uniformly in f G (0,1), and for every x G M, which proves that ||/||c‘^ 

Ia"^- 


< 


Step 2. ^ A®". Let / G be given. Using the decomposition of the identity 

provided by Calderon reproducing formula 


f = e-^f+l 


/' 




dt 


we first deduce that / is bounded, with 


f 


dt 


00 < ll/llc- ( 1 + / j] ^ WJWc- 


Moreover, for any two points x,y, with 0 < d{x,y) < 1, we have 

"1 . . ... s . ^ ' dt 

'0 ^ ' '' / J t 


fix) - f{y) = {{e ^f){x)-{e '^/)(i/)}+ {(QrV)(a:) - (gfV)(y)} j 

= [e~^f{x) - e"^/(y)| + | ix) - 


+ / {(qS'V)(x)-(q 


( 2 ) 


dt 

t 


















17 


One can use the Lipschitz regularity (Lip) of the heat kernel to bound the first 
term in the above sum, giving 

- e~^f{y) < j \pi{x,z) - pi{y,z)\\f{z)\y,{dz) 

^ d{x,y) 


As similar bounds hold for Le ^, by analyticity of the heat kernel, the second term 

admits a similar upper bound. Let now focus on the third term, using a similar 

reasoning and noting that So, for d{x,y) < y/i, we can write 

2 2 

(qPV) (a:) - (y) < f Kq(2){x,z) - K^i2){y,z) K^(2)f{z) p{dz) 

J -n TT TT 


< 


d{x,y) 

Vi 




< d{x,y) I 


Vi 

If \/t < d{x,y), then we directly have 


(QfV)(^)-(QfV)(y)|<||QfV 


< t 2 


Hence, 




< 


a dt V 
/ ^'T+/ 

'0 y Jd{x,y)^ 


d{x,y) 

Vi 


t 2 . 


dt 


< d{x,yY wjiico 
since a S (0,1). Consequently, we have obtained 

\fix) - f{y)\ < d{x,yy\\f\\c- 

uniformly for every x Y U with d{x,y) < 1, so indeed ||/||a'^ 


< 

('N_/ 


o 


Our main example of a distribution with negative Holder exponent a will be given 
by typical realizations of a (possibly weighted) noise over {M,p) - see Proposition 15.41 
To prove that regularity property, it will be convenient to assume that the metric 
measure space (M, d, p) has the following property, called Ahlfors regularity. There 
exists a positive constant ci such that 

V{x, 1) > Cl, 

for all X G M, which, by the doubling property, implies that we have 
(2.10) V(x, r) > cir'^, 

for some positive exponent u, all x G M and all 0 < r < 1. The constant is d on 
a d-dimensional manifold equiped with a smooth measure. This is a relatively weak 
assumption that is essentially satished in a Riemannian setting for closed manifolds 
without boundary and injectivity radius bounded below by a positive constant or in a 
sub-domain of the Euclidean space provided that the boundary is Lipschitz. Under that 
additional non-degeneracy assumption on the volume measure, we have the following 
Besov embedding, proved in Appendix lAl 


























Lemma 2.6 (Besov embedding). Given —oo < a < 2, and 1 < p < oo, we have the 
following continuous embeddings. 


B. 


p,p 


B. 


p,oo 


(J - 

D P 


= C° 


Besov embedding can be used in a very efficient way to investigate the regularity 
properties of random Gaussian fields, as will be illustrated in section [5l3l 

Remark. Let us point out here that our Holder spaces , with a < 0, coincide in 
the Euclidean setting with those used by Hairer [35], and, roughly speaking, defined 
in terms of scaling properties of the pairing of a distribution with a one parameter 
family of rescaled functions. Indeed, on the Euclidean space it is known that to define 
Besov spaces or Holder spaces through Littlewood-Paley functionals, we may chose any 
good Fourier multipliers satisfying suitable conditions; the latter are satisfied by the 
derivatives of the heat semigroup. So our spaces correspond to the standard 

inhomogeneous spaces defined by any Littlewood-Paley functionals. From wavelet or 
frame characterization [52], we can conclude that our Holder space coincides with those 
used in [35] or [38] , 

Before turning to the definition of our paraproduct, we close this section with two 
continuity properties involving the Holder spaces , which we shall use in the sequel. 


Proposition 2.7. For any a G (—oo,2), and every integer a > 0, we have 


p(a)f 


< 

OO 




Proof - We have by construction P^^"^ = (l + aiL + • • • + Ua-iL'^ ^)e for some 
coefficients ai,..., aa-i. As we have by definition ||e“^/||^ < ll/llc'^) and L^e~^ = 

, for .£ = 1... (a — 1), we see that ||L^e“'^/||^ < ||/||c<^, since we have seen 
above that we can choose the parameter a in the dehnition of the Holder space. 

[> 


Proposition 2.8. For a G (—oo,l), we have 


sup t 2 

* 6 ( 0 , 1 ] 


(Vtr)(e-*V) 


< 


c<^- 


The same conclusion holds with any of the operators in the role of e 
Proof — Given t G (0,1], use Galderon reproducing formula to write 

(Vtr)(e-*^/) < [' (Vtr)(e-*^QW/) -+ (Vtr)(e-('+‘)^/) 

Jo 'S 

We divide the integration interval in the above-right hand side into (0, t) and [t, 1] 
to bound that term. For s < t, we have e“2^ ^1+i’ 

Proposition 12.31 to get 


-tL 


(Vtr)(e-*^QW/) 


< _ 

OO Z 




Q'-H.f 
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Similarly for t < s, then e = e 2 Q^s]_., and we have 

2 2~*~^ 


(Vtr)(e-^^gW/) 


< 


1 

t \ 2 


{,/ir)(.-l'^) 


<3-+</ 


< - 


t\ 2 


S2 


Similar computations give the estimate 

(Vtr)(e-('+‘)^/) < Vtll/llc- 

cx> 

We conclude by integrating with respect to s € (0,1), using here the fact that 

CJ < 1. 

O 


3 _ 

Paraproduct and commutator estimates in Holder spaces 

3.1 Paraproducts based on the semigroup Bony’s paraproduct machinery has its roots 

in the Littlewood-Paley decomposition of any 
distribution / as a sum of smooth functions A,/ localized in the frequency space, so a 
product fg of any two distributions can formally be decomposed as 

(3.1) fg = Y,^^f^l9= Y. Y A./A,5=:(l)+n(/,5) 

i,j \i-j\>2 

into a sum of products of two functions oscillating on different scales, and an a priori 
resonant term Il{f,g). This decomposition draws its usefulness from some relatively 
elementary a priori estimates that show that the term (1) above makes sense and is 
well-controlled under extremely general conditions, while the resonant term n(/, g) can 
be shown to define a continuous map from x to provided a + f3 > 0. These 

estimates rely crucially on some properties inherited from the very definition of the 
Littlewood-Paley blocks as Fourier projectors. These properties cannot be grasped so 
easily in our semigroup setting; however, we shall use the operators and y/tT 

or {tL)P^°'^ as frequency projectors, with projecting on frequencies lower than 
or equal to ^“2^ and ,y/tT or {tL)pj:°^ as localizing at frequencies of order t“2. 
Indeed, in the torus, and working with the Euclidean Laplacian, the operator has 
for instance a Fourier transform equal to 

it is essentially localized in an annulus | A| ~ t“2. Similar explicit Fourier pictures for the 
other operators can be given in the setting of the torus. This ’frequency’ interpretation 
of these operators will be our main guide in the definition of our paraproduct given 
below. This paraproduct will depend on a choice of a positive integer-valued parameter 
h that can be tuned on demand in any given problem; the bigger h is, the more we can 
do some ’integration by parts’ - it will be fixed at some point, and the reader is invited 
not to bother about its value. To clarify notations, we shall repeatedly use below the 
notation / • g for the usual product of two functions. 
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Rather than starting with Bony’s decomposition (|3.ip . we take as a starting point 
Calderon’s reproducing formula 



+ A_i(/,5r), 


where 

A_i(/, 5) :=Pr 

stands for the low-frequency part of the product of / and g, and where we implicitly 
make the necessary assumptions on / and g for the above formula to make sense. 

Guided by the above heuristic argument about the role of the operators ,Q'f\ 
etc. as frequency projectors, we decompose the term involving the product of 
and g, by using the definition of the carre du champ operator F 
L{4>1 ■ h) = L{(j)i) (j)2 + L{4>2) 4>i - 2T{(j)i,(j)2) 


and write 


Qib) 

= (Pf V • (tL)Pf ^5) 

=:Bg{f) + Bf{g) + R{f,g). 




If one rewrites identity (|3.2p under the form 

Hf 

f9=- {(1) + (2) + (3)} Y + 

with obvious notations, this suggest to decompose it as 

f9 = J^ ({(1) + ^gif)} + {(2) + ^ 1 ( 9 )} + y + ^-i{f,9) 

and to identify the integral of the terms into brackets in the above formula as para- 
products, and by defining the resonant term as the integral of R{f,g). This is what 
was done in m where this notion of paraproduct, introduced in [ 10 ] . was shown to 
have nice continuity properties in Holder spaces C", provided one deals with positive 
exponents a. Given our needs to deal with negative exponents, a refinement of this 
decomposition seems to be needed to get some continuity properties for negative expo¬ 
nent as well. We thus use the carre du champ formula in each term (1) and (2), and 
decompose (1) under the form 

(tL)pf) (qS'"'V • Pt^'^9) + {2Pf ^ (tr) (gf-'V, Pt^'^9) - Pt^'^ (gf "'h • {tL)P,^^^g) } 

=-.Ag{f) + S{f,g), 


with S{f,g) the sum of the two terms into bracket, and 


{2)=Af{g) + S{gJ). 

Note that the functions Af{g),S{f,g),... all depend implicitly on time. This decom¬ 
position leads to the following definition. 


21 


Definition. Given an integer b >2 and / G Use(oi)^^ 9 ^ define their 

paraproduct by the formula 




= 1 ^ (gf-'V • p^^g) + • pf 5)} 


(if 

t 


The well-defined character of this integral is proved in proposition 13.21 below. With 
this notation, Calderon’s formula becomes 

fg = nf if) + nf (g) + n('’) (/, g) + A_1 (/, g) 
with the ’low-frequency part’ 


and the ’resonant term 

fi 


dt 

t 


n(') (/, g) = ^J^ g) + ^(5, /) + Pif, g)} 

= {-p^^ (gf"'V-(iP)pf^5) +2pf^r(VfgS'"'V,Vfpf^9)} j 


+ 


-f 

lb Jo 

1 


-K 


(b) 




- w / ^g* 

76 Jo 




{mPlf-Qf-'lg) 
(ViPlf. sTtPlg) J. 


+ 2Pf 

dt 




(if 

T 


Note that we have = Id — A _i(-, 1), as a consequence of our choice or renormal¬ 

izing constant. 

3.2 Paraproduct estimates We prove in this paragraph the basic continuity estimates 

satisfied by the maps defined by the low frequency part, the 
paraproduct and the resonant term. 

► Proposition 3.1. Fix an integer b > 2. For any real numbers a,l3 and every positive 7, 
we have 

(3-3) ||A_i(/,g)||^^ < ||/||c“||ff|lc/3- 

for every f £ and g G . 

(gi“0 

V /0<S<1 


Proof — Consider the collection 


for a large enough integer o > 7. Then 




Since s < 1, we have gi“^p/^^ = s“e with the operator bounded 

on L°°. We obtain the conclusion from Proposition 12.71 as we have 

||q7a-i(/,<,)L < .ip<*>/Li|F79L 


<s"ll/lbll9llc»- 


O 
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The continuity properties of the paraproduct are given by the following statement; 
they are the exact analogue of their classical counterpart, based on Littlewood-Paley 
decomposition, as can be found for instance in the textbook [8] of Bahouri, Chemin 
and Danchin. 


Proposition 3.2. Fix an integer b>2. For any a G (—2,1) and / G we have 
• for every g G L°° 


(3.4) 


(3.5) 


n®(/) 


< 


ll5l|oc||/lk 


for every g £ with /3 < 0 and a + /3 G (—2,1) 


n®(/) 


< 


Qa + P 


Il9llc/3||/Ilc“- 


The range (—2,1) for the regularity exponent can appear as unusual since in the 
standard Euclidean theory such continuity properties hold for every a G M. However, as 
explained in section ITTI the restriction a < 1 comes from our optimal / minimal setting 
where we only assume a gradient estimate on the heat kernel. The restriction a > —2 
can be explained as follows. In the Euclidean theory, nice Eourier multipliers can be 
used to have a ’perfect’ frequency decomposition and the study of paraproducts mainly 
relies on the following rule: the spectrum of the product of two functions is included 
into the sum of the two speetrums; this comes from the group structure through the 
Eourier representation of the convolution. In our setting, our frequency decomposition 
involving the heat semigroup is not so perfect and the previous rule on the spectrum 
does not hold, at least not in such a ’perfect’ sense. That is why the new limitation 
a > —2 appears; it is inherent to the semigroup approach developed here. No such 
limitation holds in the more restricted setting developed in [7]. 


Proof - Recall that 
1 


nf (/) = j 

Given s G (0,1], consider ^^I[g{f). Eor s < t, we use that 
= (I)'"' (tL)'pf and 


-V„-sL 


and for t < s that 

gi'-i)(tL)pf = and 

Hence, with the uniform L°°-boundedness of Qt, Pt operators, we have 


g('-i)nW(/) 


+ 


/■ 


< / - 

oo Jo s 






t 

H— 

oo S 


{tL)pj;^'^f 


pF’^g 


dt 

oo t 


b-1 


gf-'V 


P^^g 


+ 


3-1 


(tL)pfV 


P^^g 


dt 

oo t 


Since / G we have 


gf-'V 


Moreover, if 5 G L°° then 


+ 


(tL)P®/ 


< 12 


C“- 


P?g 


< 
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and if 5 G with /3 < 0 then 




< / 

°o Jt 




du 


+ 


CXD U 


P^ig) 


< 




\\g\y ^ p lbllc/5- 


We deduce the following bounds as a consequence. 
• If G L°° then 




< 


^ — 

< .S 2 


dt 


J ‘W + , 


^ /S\C a dt 

- t^ — 


t 


since a G (—2,1) and c > 1. This holds for every s > 0 which yields ()3.4p . 
If G with a + /3 G (—2,1) then 




< 

rs_/ 


01+13 

< S 2 

rs_/ 


t \ 0+13 dt 


,, . . 

^ " y ) ll/llc“lbllc/5 


since 2(6 — 1)>1>Q! + /3> —2. This holds for every s > 0, which yields 
(1331). 

[> 

Proposition 3.3. Fix an integer 6 > 2. For any a, /3 G (—oo, 1) with a + /3 > 0, for every 
/ gC“ and 5 gC^, we have the continuity estimate 




< 




Il/lbll5llc/^- 


Proof - We recall that 




"1 


7b io 


-y 


(b) 


(gf-'V • ^9) + (Vtgf "'V, Vip^^g 


+ 


+ 




+ 2P®r 


dt 

t 

dt 




lb Jo 


dt 

t 


Consider the function ^'^Il^^\f,g), for every s G (0,1]. It is given by an integral 
over (0,1), which we split into (/) an integral over (0, s), and {II) an integral over 
(s, 1). Since / G C", the use of Proposition 12.81 with a < 1, yields the estimate 




+ 


ytr(gf-'V) + + Vtr(gf-'V) 


(b) 


lib-1) 


+ 


Vtr(pfV) 


< t2 


a similar estimate holds with g in place of /, and /3 in place of a. Using the uniform 
L°°-boundedness of the different approximation operators, we get for the first part 




WfWMc^ 


0+13 
< s 2 


z°‘\\g\\cP 
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where we used the strict inequality a + /3 > 0. For the second part, we observe 
that for t > s then 

and 

So we get for the second part 

^ a + /3 /s \ b-i 


< 3 f-‘'(n)L S (/ (i) 


y lll/IbNIc/^ 


< s""''ll/llc“||5'llc/9> 


using the fact that 2(6 — 1) > 2 > a + /?. 


O 


3.3 


Commutator estimates Recall the discussion of the paracontrolled calculus approach 

to the study of the parabolic Anderson model equation given 
in the introductory section 11.21 We have at that point the paraproduct and diagonal 
operators in hands; we deal in this section with the fundamental commutator estimate 
introduced in [32]. Readers familiar with the basics of stochastic analysis will notice 
the similarity of this continuity result and the rule satisfied by the bracket operator in 
Ito’s theory 



Pd{N,py, 


this is obviously not a coincidence. 


► Proposition 3.4. Consider the a priori unbounded trilinear operator 

C{f,g,h) :=nW(nf (/),h) -gn('')(/,/r), 

on Sy Let a,/ 3,7 be Holder regularity exponents with a G {—1,1), (3 G (0,1) and 7 G 
(— 00 , 1 ]. If 

0 < a + /3 + 7 and a + 7 < 0 
then, setting 6 := (a + /3) A 1 + 7 , we have 

(3.6) ||C'(/,5,^)||c^ ^ ll/llc“ hWcP Iloilo, 

for every f £ ,g £ and h £ C'^] so the commutator defines a continuous trilinear 
map from x x to . 

Proof - Note first that the paraproduct 13g{f) is given, up to a multiplicative constant, 
by the sum of two terms of the form 

A{f,g) = 1^' Q][QU-Vtg)j, 

and the resonant part n(/, g) by the sum of five terms of the following forms 

(3.7) n{f,g) = {VtV^f, VtVh) j, 

or 

r ^ (jf 

n{f,g)= / Vt{{tL)Vtf-Qtg) 

JO t 

or 

/■i gf 

nf,g)= / rHQtf ■ {tL)vh) 

Jo ^ 


where the operators 
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• Qt, Qi s-re of the form (tL)^ ^p{tL) e with a polynomial function p, 

• Vt^Vl are of the form p{tL) with a polynomial function p. 

Note also that terms of the for 'ip{tL) are a posteriori of the form (j)(tL). So it 
suffices to focus on a generic term of the form 

D{f,g,h) :=TZ{A{f,g),h) - gTZ{f,h) 

and prove the continuity estimate (|3.6p for it. We focus on the case where TZ has 
form ()3.7p . the treatment of the other cases being similar and somewhat easier. We 
split the proof of the commutator estimate ()3.6h for D in two steps, and introduce 
an intermediate quantity 

S{f,g,h) := iyi{r{VtV^f,ViV^h) • Ptg} j 

for which we shall prove that we have both 

(3.8) \\gn{f,h)-S{f,g,h)\\^, < \\f\\c<^ \\g\y ||h|b 
and 

(3.9) \\D{f,g,h)-S{f,g,h)\\^s < ||/lb \\h\h ■ 

Step 1 — proof of (j3.8p . We first prove a weaker version of the continuity estimate 
(13.81) . under the form of the inequality 

(3.10) \\gn{f,h) -S{f,g,h)\\^ < ||/||c« \\g\y \\h\y. 

As a start, remark that we have 

(3.11) {gn{f,h) - S{f,g,hfj{x) = j Vl\T{VtV^ f .VtV'fh) ■ {g{x) -Pt5 r)|(x)y, 

for /r-almost every x G M. Since g G C^, with /3 G (0,1), we have 

\\Ptg-gU< / IIQ.slloo- 
Jo s 

^ (i) ~ 

so we have 

\Ptg{y) - g{x)\ < \Ptg{v) - g{y)\ + \g{y) - 5(a^)| 

^ d \\g\y +d{x,yf \\g\\cP 

< +d{x,yf'^ \\g\y, 

for every x,y £ M. Coming back to equation (j3.11ll and using Gaussian pointwise 
estimates for the kernel of V^, together with Proposition 12.81 and (|2.6p . we have 

pi (T{ViV?f, Vtvfh) . {g{x) - Ptg)){x) 

no greater than 

{Im V{x Vi) Il^^(^"^)llooll^r(p3/i)||^ 

^{Im V{x Vi) +d{x,y)f^'j dp{y)^ ft ii ||/||c“||£(||c/3||^llo 
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The continuity estimate (j.S.lOh comes from integrating with respect to time, taking 
into account the fact that a + /3 + 7 > 0. 

Let then estimate the regularity of gTl{f,h) — S{f,g,h). For G M, with 
d{x, y) < 1, write 


(^ 9 ix)'Jl{f,h) - S{f,g,h)^{x) - {g{y)n{f,h) -S{f,g, 
with U defined by the formula 

pd{x,y)'^ '' 


(y) =:U + V 


V. 


I (T{VtV^f, Vivfh) • (gix) - Ptg)) (x) 


-Vl T 


'{Vtvff, Vtvfh) • {g{y) - Ptg)^{y) \ j, 


and V is defined by the formula 
rl 


V. 


Jd{x,yf 


I (T{VivU, Vivfh) • {g{x) - Ptg)) (a 


-V, 


I Vivfh) • {g{y) - Ptg)){y) \ j. 


By repeating, the argument used in the proof of ()3.10p . we easily bound U by the 
quantity 

/ rd{x,y)^ c,+0+-y dt\ 

^ " jj \\f\\c°‘\\9\y\Mc^ 

<d{x,y)^\\f\\c^\\g\y\\h\\ci. 

For the second part, we use the inequality 

|B| <A + B, 

with A equal to 


P, 


'd{x,y)^ 


and 


I (F(VtnV, VtPfh) • {g{x) - Ptg)){x) 


- P 


I {T{VtPU. VtP^h) • {g{x) - Ptg)) (y) 


dt 

t 


B:= \g{x)-g{y )\• Pl(T{VtPU,VtP^th)){y) 

Jd{x,y)'^ ^ 

The last quantity is bounded by 
B ^dix,yf \\g\y [ 

Jd{x,y)^ 


< d{x, yf ll/llc- \\g\\cP\\h\\c^ I y 

^2 t 


Vir{v^f) Vir{vfh) 

00 

r 

d{x,y)' 


dt 

t 

dt 

00 t 


< 

rsj 


<!(*,!;)* ll/ll<;.|l9llco|l'>lb- 
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For the quantity A, we use the Lipschitz regularity (Lip) of the heat kernel to get 
the upper bound, 


d{x,y) 

d{x,y)^ \JM '\/tF(x, "x/t) 


exp —c 


d{x, zY 


\g{x) - Ptg{z)\y{dz] 


VtT{V^f) VtT{V^h) 


dt 

T 


< 


d{x,y) 


^dix,z)^ 


(d{x,zf fj.{dz)t’^ 2 '' I ||/|lc-||5llcH|h||c 


d{x,y)^ JM y/t\^ y/Tj 

d{x,y) dt 

-- 1 2 - 

cZ(x,y)2 y/t t ^ 

<d(x,yr+^+^ii/iic.iiffiic.ii/^iic., 


< 


ll/llc<»||ff||c.||/^ll 


C'l' 


where we have used the fact that a + j5 + ^ G (Ojl)- The combination of all the 
previous estimates yields 


{gT^{f, h)-S{f,g,h)^{x) - \gn{f, h) - S{f,g,h 


<\U\ + A + B< d{x,y) ||/||C“||5llc/3|l'^llo, 
which concludes the proof of the continuity estimate (|3.8I) . 

Step 2 — proof of (13.91) . Given the collection [Qr ■= of operators, 

V /re(0,l] 

we need to prove that we have 


(3.12) 


r* 2 


Qr(TZ{A{f,g),h) -S{f,g,h)^ 
for every r € (0,1], and where TZ[A{f, g), h) — S{f,g, h) is equal to 

( 3 . 13 ) (^Vtyy^QliQlf-Pig) ^-p.g.p^f^ ,VtPfh}j J. 

The notation may be confusing and we need to be careful: when F acts on the 
product Ptg ■ Ptfj h is thought to act only on the variable of V^f, with the 
variable of Ptg frozen. We shall bound above the absolute value of the F term in 
the integral, which is of the form F(p, g), by F(p)F(g) - recall we write F(p) for 
y^T{p,p). Set for that purpose 

At{f,g) := ViT ^ V^Q\ {Q^f • P^g) ^ - PtgP ^. 

We have for almost every x G M 

Mf,9){x) < VtTp'^ y Q\ {Qlf • P^g) (x) ^ - Ptg(x) ■ (x) 

< r VtrP^Ql [alfiP^g - Pt9ix))){x) ^ + \Ptg{x)\Vtr[P^Pif]ix), 


where we used the property 


L 




= Id-iPi, 
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for some Vi operator. As in Step 1, the fact that /3 is positive, which implies the 
inequality 

vlgiy) - Ptg{x) < V^giv) - giv) + laiv) - g{x)\ + \g{x) - Ptg{x)\ 

< (^s2 +t2 +d{x,yf'^ WgWcP < (max{s,t)2 + d{x,y)^'^ hWci^- 

Moreover, it follows from Lemma lA.51 in Appendix A ~ about the composition of 
Gaussian pointwise estimates, that the operator y/tT has pointwise Gauss¬ 
ian estimates at the scale maxfs, t) with an extra factor ( j ^ ■ gQ Qj^g gg^g 

’ '' \ max(s,t) J ’ 

r := max(s,t), we have 
VirpfQl (Qlf{Plg - Pt9(x)))(x) 

1 „ d(x,y)^ / /3 


< 


< 


< 


/ min(s,t) 

\max{s,t) J \JmV{x,Vx)' 
/ min(s, t) 

\max(s, t) 

f min(s,t) 


+d{x,yf'j g{dy)> \\Qlf\\^\\9\\ei^ 


||/||c“||ff||c/3 


max(s,t)2 s5 ||/||c-||3||c/9- 


\max(s, t) ^ 

Integrating in s, and taking into account the fact that a > — 1 and a + P < 1, we 
obtain for At{f,g) the estimate 

rt 






2 .P a ds , f t\‘^ P a ds /-[ . 

S2S2 _ + k ll/llc^lbllc/? 


^+0 


2 ll/llc“lbllc/5- 

Observe that in the case where a + /3 > 1, we get 

\\Mf,9)\\^ < t^WfWc^-hWcp^ 

Goming back to identity (13.131) . we have 

7^(M(/, g),h)- 5(/, g,h)\< Vj (^Mf,g) • VtT{Vth)) j 

and since q; + /3 + 7>0, it follows that 

P{-^if,g),h) - S{f,g,h) 


< 

rs_/ 


a + P+~i dt 
/o t 


\\f\\McP\\h\\ 


C-y 


^ ll/llc“||fi'llc/5|l^llcT. 

Moreover, taking into account that we have Q^r^V\ = jPrQt^'^ for t > r, and 
q; + /3 + 7<1, we see that the estimate (|3.1‘2I) holds true 

Ql}'>(n{A{f,g),h)-S{f,g,h))\\^ 

r\\Mf,g)\\^ Vtr{vfh) r!-\\Mf,g)\\^ Vtr{vfh) 

J Q ^ Jt ^ 

/ 


< 

rs_/ 


dt 

oo t 


< 

r\^ 


a + 0+y dt a + /3 + 7-2 dt 

t2- I rt 2 — 

IQ t 


< 




rs_/ 
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O 


3.4 Paralinearization and composition estimates Two ingredients are needed to turn 

the machinery of paraproducts into an 
efficient tool. To understand how nonlinear functions act on Holder functions with 
0 < a < 1, and to understand how one can compose two paraproducts. The first point 
is the object of the following analogue of Bony’s classical result on paralinearization 
m, while the second point is dealt with by theorem 13.61 below. 


► Theorem 3.5. Let fix an integer b>2,a£ (0,1), and consider a nonlinearity F G C^. 
Then for every f G C°^, we have F(/) G C" and 

Rp{f) 


More precisely 


F(/)-n®^)(/) 



^ 11^11(73 (^1 + ll/llc“) • 


If F G then the remainder term Rpif) is locally Lipschitz with respect to /, in so far 
far as we have 


Rpif) — RF{g)\\(. 2 a ^ I|F||c'4(i + ||/||c“ + lbllc“) 11/ — 


Proof — First using the Leibniz rule for the operator L, we know that for h G C°‘ then 

L(F(/i)) = F'{h)L{h) + F”{h)T{hf. 

Now, since the semigroup is continuous at t = 0, we have 

F(/) = limPf)F(Ff/), 


so we can write 

F(/) = - )F(pf V) dt + PWF(P('V) 

= (F(pf V)) + Ff) (gf / • F'(pf V))} y + ^F(p('V). 

Using the relation together with the chain rule 

L(F(pf V)) = F'(pf V + F"(pf V)r(Pf V)', 

we get 

gW(F(pW/)) = QS'-F(^(iL)pW/.F'(pW)) +Qf-i)(F"(pfV) •tr(P®/)'). 

Note here the identity 

pf) (gfV . F'(pf V)) =(tF)pf) (qS'-'V • F'(pf V)) - ^ • tLF'(pf)/)) 

- F'(pfV)). 


F(/) - =: (a) + (6) + (c) + {d) + (e) + (/) 


So we have 
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with 


(а) :=Pf)F(P®/), 

(б) := (F'(PfV) -^’fHF'(/)))} 7 

(c) := (F'(PfV) -^f^(F'(/)))} y 

(d) := y ^'Qf-'){F"(pfV) •tr(pf)/)'} j 

(e) := -y Pf^{QS'"'V • tLF'(pfV)} j 

(/) := -y ^'pf^tF(Qf-'V,F'(pfV)) j- 

We are now going to control each of these terms in the Holder space C^". 

Step 1 — term (a). Since / G C", we know that £ F°°, so F(P^^^^/) is also 

bounded. From Proposition 12.71 we get 


a)llc2c. < 


F(PfV) 


< 


ll/llc- 


Step 2 — terms (6), (c). The following quantity appears in these two terms 

F'(pW/) _pW(F'(/))| < ||f'(/) -F'(P®/) 

< ||f"|| 


< IIF" 


< ||f" 


P'(/)-Pf)(F'(/)) 
/-P®/| + F'(/)-P®(F'(/)) 

I oo 

-+ f <3f>(F'(/)) 

■S Jo 


I 


oo 

ds 


OO S 


- I ll/llc" + 


ds 


y 1 IlF' 


C“ 


(3.14) 


<td|F"Lll/llc"! 


we used along the way the characterization of Holder space, for Oa < 1, given by 
Proposition 12.51 to see that 

iiF'(/)iu s r'L ii/ik". 

Using this estimate (I3.14p . we deduce the following bound. Uniformly for every 
s G (0,1), we have 


Qi^hb) <( [" 

°o \Jo 


t 2 


oo—1-00 


< 


fie* 

0 t J s b t 


F'(pf/)-Ff)(F'(/)) y) ll/llcii 

oo Z y 

"ILIl/llc“ 


£«1|F"ILII/II?". 

where we used that 5- (7t) • Fhat yields 

llwlU s IF"Lil/il?., 

and a similar inequality holds also for the third term (c). 
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Step 3 — terms (d), (e) and (/). We quickly sketch the boundedness of each of 
these three terms. Using Proposition 12.81 we get a bound uniform in s G (0,1), of 
the form 




< 


< 




- 1 ) 




2 dt 


— F 
t I " 


Jo t J, t t J '' 


2 




2 

C“- 


Similarly 

ll<5i‘’(/)L;S 




dt 

oo t 


< 


/■%«*+ 

Jo t J, t t 




where we used the Leibniz rule 


^r(F'(pfV)) 


< llF"| 


Vtr((p®/)) 


< t2 IlF" 


C“- 


For the remaining last term (e), we can still nsing the Leibniz rule and get 


tLF'(pf/) < (| 

OO \ 


which then yields 


+ F" 


) ( (tL)pf / + V^r(pfV) ' ) 

/ y oo oo J 


ll«L(e)IU< /‘ 

Jo 


p. 


(fe) 




00^00 
/“I 




tLF'(pfV) 


dt 

oo t 


||F||c.|l/|l?. 

<*“||F||c. ll/ll?.. 

By combining the previous estimates, we conclude that we have 

||(<!)|lc». + l|(e)|lc>. + ll(/)IU S llfllc- ll/llc (1 + ll/llc"). 

which ends the proof of the estimate of the remainder. The Lipschitz regularity 
of the remainder term is proved by very similar arguments which we leave to the 
reader. 

t> 

Let us now examine the composition of two paraproducts. Note that for tt G and 
u G with a G (0,1), /3 G (0, a], we have uv £ . 

Theorem 3.6. Fix an integer b > 2, a £ (0,1), /3 G (0,a] and consider u £ C°^ and 
V £ . Then for every / G we have 


ni')(nW(/))-nS(/)GC“+^ 


with 


nW(nW(/))-n(tH/) 


< 




C“ lkllc“ ll^^llc/9- 
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Proof - We leave a detailed proof to the reader and we just sketch it, since it is similar 
and easier than the proof of Theorem 13.51 Following Proposition 13.21 we know 

that the two terms and nl^2(/) belong to The idea is to use the 

C^-regularity of v to gain the same regularity in the difference. 

Indeed, adopting the notations used above, the paraproduct Ilg^^ (/) is given, up to 
a multiplicative constant, by two terms with the form 


Hf,g) = I' 


where in Qj and Qf we have at least a term (tL) to the power 1. Let us focus on 
this form. Then we have 


/(/(/, n), n) = Ql [qIqI {Q^f • • V^u 


ds dt 
st 


and 


l{f,vu) = £ Qj(Q^Ql {Qlf)-rhnv)) 


ds dt 
st ’ 


where we have used the normalization = Id, which means here that /(/, vu) = 
l),nu). Then using the C^-regularity of v and the fact that Ql involves at 
least a power 1 of (tL), one can check that uniformly in s,t £ (0,1)^, we have 


qUq^qI {Qlf • rlv) • Piu) - aUajal {alf) • vl 


[uv 


oo 


max(s, t) 

So integrating in s £ (0,1) yields for a + /3 < 2 


, minis,t) a p . 

~ TTTTTTT®" (s + t) 2 ||/||c“ lkllc“ • 


I 


' "qUq'^QI {Q‘if-'Plg)-Vlu) - Qj(Q^tQl {Q^f) ■vHuv)'^ 


OO S 


a: + /3 


2 ll/llc“ll^^llc/^lkllc“- 

Then as previously we check that for every r £ (0,1) we have 


< 

rs_/ 


min(T,t) ^(^+ 0)/2 dt 


a+P 

< 7-2 


Q max(T, t) t 

c<-\\v\\cp\\u\\c^ 


C“ ItIIc^ ||'^||c“ 


since a + /3 < 2. That allows us to conclude that 




< 


C°-+P 


ll/llc“lbllc/3|l«l 


|C“ • 


o 
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3.5 Schauder estimates Proposition 13.81 gives an elementary proof in our setting of a 

Schauder-type estimate about the regularizing character of the 
convolution operation with the operators Its paracontrolled analogue, given in 

section W?2\ provides a crucial ingredient in the study of parabolic singular PDEs, from 
the point of view of paracontolled distributions. 

► Definition 3.7. For a G (0,2) and T > 0, we set 

:= I / G ||/||ctC-= sup ||/(t)||^„ < oo I 

1 t&[ 0 ,T] I 


and 


||/(t) — f{s) 

/GL-([0,r]xM), ll/ll . .= sup 

' s^t t-S2 

0<s,t<T 


< oo 


We then define the space 


f ;= C7 tC“ n (P| L" 


The resolution operator TZ is formally defined by the formula 

n{v)t:= [ ds; 

Jo 

it implicitly depends on the parameter b. 

Proposition 3.8. Consider an integer 6 > 0 and a regularity exponent /3 G M. For every 
positive finite time horizon T, and every v G CtC^ , then TZ{v) belongs to and 

\\T^{v)t\\c !^+2 < (^ + T) sup ||i^(s)||c/3, 

se[o,d 

for every t G [0, T]. Moreover if —2 < /? < 0 then we also have 


|P(u)|| 


/3+2 

L° 


< 


’WctCP- 


Proof — We consider another integer c > -^ + 1 and a parameter r G (0,1]. Then 


Q^\n{v)t) = pelvis) ds. 

Jo 


We have 


(3.15) 




< 


T + t — S 


(t + t- s)^\\v{s)\\cp. 


Indeed, if t — s < r then we only use that Qr and Pf_g commute with the L' 


boundedness of Pt^g to have 


Qi^'>Pt\v{s) 


< 




V[s, 


< 


T2 ||u(s)||c/3. 


In the other hand, if t — s > r then we use that 
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and we conclude similarly with the L°°-boundedness of e and the property that 


{{t - s)ypi’^\v{s) 


< / 

oo Jt- 


< 

rsj 


< 

rsj 


{{t - s)yQi’’K{s) 


dr 
oo r 



' t — S 



t — S 


' t—s 




lk(s)llc/5 


dr 
oo r 
dr 
r 


That concludes the proof of (|3.15l) . So by integrating, it comes 




T + t — S 


{t + t- s)^ ds \ sup ||i;(s)||c/3 
se[o,t] 


<r 2 +i sup ||w(s)||^^. 
se[o,t] 

This holds uniformly in r S (0,1] and so one concludes the proof of the first 
statement with the global inequality 

l|KW.L< 


(■ ■ ■ )^'S f lbllctC/5 ^ ^IbllcrC/S- 
we have 


For the second statement, we note that for s 

n{v)t - n{v)s = (pE - Id) {n{v)s) + pyi{v{r)) dr 

= — [ Qy'>'Jl{v)s — + [ PyX{v{r))dr. 

la JO f Js 


We have 


llTew.lU 


< 

r\j 


and since /3 < 0, we also have 

^yib) 
f s 


Pt-r{v{r)) dr 



pE E(^)) 
pE E(^)) 


dr 


dr 


where we used ^ + 1 G (0,1). 

Corollary 3.9. For a fixed integer 6 > 0 and a G (0,2), we gave 

(^)||£:(“) ~ (I + ^) II^IIctC“-2) 

[0,T] 


\> 


uniformly in T > 0. 
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Remark 3.10. Observe that in Proposition \3^ the weight (1 + T) can he weakened, 
up to a little loss on the regularity exponent. Indeed, the exact same proof allows us to 
prove 

||7^(^;)t||^^+2-2e ^ sup ||^;(s)||c/3 
se[o,4] 

and 

for any e € ( 0 , 1 ); so we have 

||^(^)L(-) ^ ll^llc'TC“-2+2e• 

[0.T] 

We refer the reader to Provosition \5.^ for a detailed proof of a more difficult statement, 
where we show how we can improve the bound (l+T) up to a small loss on the regularity. 


Paracontrolled calculus 

The idea of paracontrolled calculus, such as introduced in [32], has its roots in Gu- 
binelli’s notion of controlled path m- The latter provides an alternative formulation 
of Lyons’ rough paths theory [sniEi] that offers a simple approach to the core of the 
theory, while rephrasing it in a very useful Banach setting. Let us have here a glimpse 
at this field, as a guide for what we shall be doing in this section and the next one. 
We refer the reader to [21| for a very nice and pedagogical introduction to the subject, 
assuming only here that she/he knows only the very definition of a (weak geometric) 
a-Holder rough path, for some g < a < |; see also [6]. Let just mention that these 
objects are nothing else than objects that play the role of the collection 

i ht — hg , / / hu®hj. dudr J 

V Js Js J 0<s<t<T 

of the increments of an M^-valued control /i in a controlled ordinary differential equation 
xt = yi{xt)h\, together with its second order iterated integral. These quantities are 
precisely what appears when making a second order Euler-Taylor expansion of the 
solution to the equation. Rough paths are such kind of objects for which h is too 
irregular to make sense of the product dhu <S> dhj.] typical irregularity is a-Holder, with 
a < like for Brownian motion. 

Assume we are given an M^-valued (weak geometric) a-Holder rough path 

X = [{Xts,^ts)) 

with Xtg € and S (8) recall that X is a-Holder and X is (2a)-Holder. 
Let also cr = (Vi,..., W) G C'^(M‘’*, L(M^, M'’*)) be given, with each column V) seen as a 
vector field on R'’*. Following Lyons, an R'^-valued path x, is said to solve the rough 
differential equation 

(4.1) dxt = u{xt)X.{df) 

if one has, for every smooth observable /, the second order Euler-Taylor expansion 

(4.2) f{xt) - f{xg) = Xlg{VJ){xg) + Xi^{V,Vkf){xg) + 0{\t- 

for all 0 < s < t < T, for some constant a > 1. Note that if Xfg = ht — hg, and 
Xtg = Jg{hr — /is) 0 dhr, for some R^-valued control h, equation (|4.2p is nothing but 
a second order Taylor expansion for the solution to the controlled differential equation 
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xt = a{xt)ht. Gubinelli’s crucial remark in [3T] was to notice that for a path x, to 
satisfy equation (US), it needs to be controlled by X in the sense that one has 

(4.3) Xt - Xs = x'^Xts + 0{\t - 

for some L(M^, M'^)-valued a-Holder path x',, here x^ = a{xs)- This set of paths has a 
natural Banach topology. The point of this remark is that, somewhat conversely, if we 
are given an L(M^, M'^)-valued a-Holder path {z,z') controlled by X, then there exists 
a unique M'^-valued path y, whose increments satisfy 

yt-ys = ZsXts +ZgXts + 0 {\t- s|“), 

for some exponent a > 1. Note that we indeed need the full rough path X to define that 
path, and not just X. With a little bit of abuse, we write Zg X(ds) for that path y, ~ 
this path depends not only on z but rather on (z,z'). Given an M'^-valued path {x,x') 
controlled by X, and cr sufficiently regular, the L(M^, ]R‘^)-valued path Zg '■= cr{xg) is 
controlled by X, with derivative z' = a'^^(xg). With Zg = a{xg) and x, satisfying the 
first order Euler-Taylor expansion (14.31) , the above second order Euler-Taylor expansion 
for y is nothing else than the right hand side of (14.2p . with the identity in the role of 
/. It can be proved that the rough integral Zg X{ds) depends continuously on (z, z') 
and X in the right topologies - this is the main selling point of rough paths theory. So, 
for a path x, to solve the rough differential equation (14.ip , it is necessary and sufficient 
that it satisfies 

xt- Xg = a(xr) X(dr), 

for all 0 < s < t < T, that is, x, is a fixed point of the continuous map 

X, I—)• / cj(xr)X(dr), 

Jo 

from the space of paths controlled by X to itself. The well-posed character of equa¬ 
tion (SH) is then shown by proving that this map is a contraction if one works on a 
sufficiently small time interval. 

Our present setting will not differ much from the above description. We aim in the 
sequel at solving equations of the form 

{dt + A)u = F{u)C, 

for some distribution C. Comparing this equation with the rough differential equation 
(BU), the role of the rough path X will be played in that setting by a rough distribution 
C= withC in the role of dXt ~ or dh, with somehow in the role of dXt 

- ox dh ® dh, and the operator (9^ -|- A) in the role of the time derivative operator 
The elementary insight that the/a solution u should behave at small space scales 
like the solution Z to the equation (dt + L)Z = is turned into the definition of a 
distribution ’’controlled by” Z, such as given below, using the paraproduct as a means 
of comparison, for writing a first order Taylor expansion of u similar to identity (|4.3p - 
compare this ansatz with the fact that a solution to a rough differential equation should 
be controlled hy X = f dX. The crucial point of this definition is that one can make 
sense of the product F(u) Q, in that controlled setting, which provides an analogue of the 
definition of the product (j{xg) X.{ds), given by the right hand side of identity (|4.2I) - see 
theorem an To run formally the above argument, we shall need to see how controlled 
distributions are transformed by a nonlinear map; this is the content of theorem 14.21 on 
paralinearisation. Some problems intrinsically linked with the multidimensional setting 
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of the problem are dealt with in section 14.21 where a version of Schauder theorem is 
proved for paracontrolled distributions. 

A last look at section o may provide a helpful guide for this section, before pro¬ 
ceeding. 


4.1 Paracontrolled distributions We fix throughout that section and the next one an 

integer b >2. 

► Definition. Let a G (—2,1) and /3 > 0 be given, together with a reference distribution 

Z G C". A pair of distributions (/, /') is said to be paracontrolled by Z if 

(/,/')“ ■.= f 

In such a case, we write {f,f') G Ca{Z) and define the norm 

||(/./')|lc^= ||(/,/')«IU« + ll/'IG. 

If Z G is a possibly different reference distribution and {h,h') G Cq(Z), we set 

d^,{{f,n,{h,h')) := Wf'-h'Wc, + ||(/,/')« - 

Note that this choice of norm allows to compare paracontrolled distributions associated 
with different model distributions X and Z. Following the terminology of [32], the 
function /' is called the derivative of /, and the term the remainder; one 

should think of the decomposition 

f = uf{X) + if,f)^ 

as a kind of hrst order Taylor formula for /, in terms of regularity properties. The 
notion of derivative depends of course on which model distribution is used. As a first 
step towards completing the above program, the following statement gives an analogue 
in our setting of the right hand side of identity ()4.2p defining (j{xs) X(ds) in the rough 
paths context. It is motivated by the following simple regularity analysis based on 
propositions 13.21 and 13.31 giving regularity conditions for the well-posed character of 
paraproducts and resonnant terms. Given / G and Q £ , with 7 < 0 <a< 1, we 

have from Calderon’s identity the formal identity 

/C = nf (C) + nf)(/) + n(')(/, c) + a_i(/, c), 

where the only term that is potentially undehned is the diagonal term n(^)(/, C). If 
however, / is controlled by Z, with derivative f G we can write 

n(')(/,C) = n(')(n(';(z),c) +n('’)((/,/0»,c), 

with ^(/, well-defined if (a -t- /3) -|- 7 > 0. So, writing 

nW (uf) (Z), c) = C(Z, f,g) + f (z, c), 

we finally see that, in the paracontrolled setting, the only undefined term in the above 
a priori decomposition of fg is the term Il^^\ZX)- The following theorem turns that 
elementary regularity analysis into a constructive recipe for dehning fff 
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► Theorem 4.1 (Product operation in a paracontrolled space). Let7<0</3<Q:<1 

be regularity exponents such that 

a + 7<0, a +/3 + 7 e (0,1). 

Let {Zn)n>o be a sequence of smooth functions converging to Z in C", and {Cn)n>o be a 
sequence of smooth functions converging to C, in C"^. Assume that Il^’^'>{Zn,Cn) converges 
in to some limit distribution write C foi" the pair (C, - call it an enhanced, 

or rough, distribution. 

• Given a paracontrolled function {f,f') G Ca{Z), the formula 

(4.4) if, f) ■ c := n(')(C) + uf\f) + nw ((/, /')», C) + ciz, f, c) + /' 

defines an element of which satisfies the estimate 

(4.5) ||(/,/') • c-n^(C)IU+. < ll(/>/')|lcg {llClb + ll^lbllClb + 
so ((/,/')-C,/) Gq(C). 

• Furthermore, this operation is locally Lipschitz in the sense that if (Z, C,C ) is 

another set of objects similar to (Z, (j, , and if {h,h') G Ca(Z), then we have 

the estimate 

ft')'?, 

(4.6) ^ _ 

< Cm {deg ((/, /'), (h, h')) + ||Z - Z||^, + lie - die. + ||c.+. } , 

where Cm is a positive constant with polynomial growth in 

M := max{||Z|lc., HCIIc, ll^llc“> lldic.’ llc<>+.> II(/>/')||cgd^')||cg}• 

Running backward the computations preceeding theorem 14.11 one sees that (/, f') ■ ( 
coincides with if / and C are both smooth - in which case one can choose f = 0. 

Proof — We examine the regularity of each terms of the defining identity (j4.4p . By 
Proposition [321 we have U^p{() G and 11^^^/) ^ Proposition 13.31 yields 

that ((/, /')ti, d G Applying Proposition 13.41 on the continuity prop¬ 

erties of the commutator operator, with a + 7 < 0, we see that C(Z,/',d G , 
with 6 = min(a + /3,1) -F 7. Since a + 7 < /3, then G Each term in 

formula (14.411 is then an element of C^, and (j4.5p holds since a < 1. The proof of 

the Lipschitz estimate (14.6p is left to the reader. 

O 


Let insist here on the fact that Z is not sufficient by itself to define a product op¬ 
eration, and that different choices of provide different definitions of the product 
operation. In another direction, using the paralinearization formula, we are able to 
study the action of a nonlinearity on paracontrolled distributions, giving us the equiv¬ 
alent of the elementary fact that, in the above classical controlled setting for rough 
differential equations, the image by some map cr of a path {x, x') controlled by some 
reference path X. 


► Theorem 4.2. Let 0 < /? < a < 1, a reference function Z G C" and (/,/') G Ca[Z) 
be given, together witha function F G 1^)- Then (F(/),F'(/)/') belongs to Ca(Z), 

and 


(F(/),F'(/)/') 


<l|F|lc»(l + ll(/./')Hcs)(l + l|Z|l?.). 
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Moreover, this operation is locally Lipschitz in the sense that we have, with the same 
notations as in theorem [4.11 a nd {h,h') e C^(Z), 

dc,((F(/),F'(/)/'),(F(/i),F'(/i)/r')) <Cm||F||^4 ((/,5), A:)) + ||Z - Z||^,) , 

where Cm is a constant with a polynomial growth in 

M := max|||Z||c«, ||(/, /')ilcg> II^0||cg} • 

Proof — Given / £ C“, we have F(/) £ C", since F is Lipschitz. We know that 
since F'(/) £ C“ and /' £ . Using the notations of Theorem 13.51 

we have 

F(/) - n® + i?F(/) 

with a remainder Rpif) G C Since / = + (/, /')**, we have 

+ n® ((/, /')«) 

with ((/,after Proposition 13.21 So 

F(/) - + c“+^. 

Using Theorem 13.61 we deduce that 

F(/)-<)^^^,(Z)£C“+^ 

which indeed shows that (F(/), F'(/)/') £ CpZ). We let the reader check the 
Lipschitz inequality for this operation. 

U 


4.2 Schauder estimates for paracontrolled The above definition of a paracontrolled dis- 
distributions tribution is adapted to a time-independent 

setting. To deal with the time-dependent 
setting needed to handle the parabolic equations considered in practical examples, we 
use an adapted notion. Recall the definition of the space given in definition 13.71 


Definition. Let a £ M and /3 > 0 be such that a + (3 £ (0)2); fix a reference distribution 
Z £ for some finite positive horizon T. A pair of distributions (/, /') £ £^ 

is said to be paracontrolled by Z if 


(«) „o(/3) 
[o,T] ^ ‘"[o,r] 


(/, /')«;= / - upiz) £ CtC“+^ n c^L°°. 


In such a case, we write (/,/') £ 


a,[ 0 ,T] 


(Z) and define the norm 


\\ifJ')\\c^ :=II/'IU + ||(/,/')" 


“a,[ 0 ,T] 




+ 11 (/,/') 


C^L° 


If Z £ is another reference distribution, and {h,h') £ Q.^ j,{Z), we set 


((/,/'), {K^)) ■■= II/' - /i'll^W + II(/,/')“ - {Kh! 


:«.[ 0 .T] 


CTC“+^nc.^ L°° 


Let us point out that the remainder (/, /')** is not a priori in p 

OL + 0 

Crp^ L°°, as the time-regularity of elements in the latter space is stronger than what 
is assumed in the above definition. 
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► Theorem 4.3. Let (5 G (0,1), a G (0, 2 — /3), and a fixed positive finite time horizon T be 
given. Given C, G , let Z be the solution on [0, T) of the equation 

CX := {dt + L)Z = C, = 0. 

Given /' G and h G , denote by / the solution to the initial value problem 

Cf = uf){0 + h, /|^^^ = /ogC“+^. 

Then (/,/') G ^ and 

||(/)/0||c^ ~ ll/ollc“+^ + (1 + r) |ll/ll£:^(i + IICIIctC“-2) + ll^lbTC“+/5-2| • 

Moreover, the map which associates {f,f') to {(, Z, f, h, fo) is locally Lipschitz. 


Proof - Since ( G CtC'^ ^ the Schauder estimates, Corollary 13.91 yield that Z G 
and / G e!^^] so we are left with checking that 

(4.7) (/, /')“ := / - UP(Z) G CtC“+^ n C^L°^. 

Let us derive an equation for this quantity 


^(/, ff =cf- cii^P{z) = n^':^(c) - cnfiz) + h 


{b). 


jib). 


= {cz) - ciE'f’{z) + h 


(b). 


c,u 


r 

{Z)^h. 


We have h G and we have seen that G and U^p(Z) G 

CtC^", so that CIiP(Z) G CtC^~‘^ ■ By studying the difference (which consist to 
commute the paraproduct with C) with introducing an intermediate time- 
space paraprodudU, such as done in |32l Lemma 5.1] - whose proof can easily be 
extended to our setting, we obtain that 


n 


c,u 


(Z) G 


where we recall that TZ is the resolution operator of heat equation. We invite the 
reader to check the Lipschitz inequality for this operation, in terms of (^,Z,f',h 
and fo- 

O 


The fact that TZ 


C,U 


iz) 


takes values in n L°° is the reason why we 


define the space of paracontrolled distributions as above rather than with the more 
natural norm Ep'P to measure the size of the remainder. 

With this result in hands, we now have all the theoretical apparatus needed to study 
some examples of singular parabolic PDEs. We have chosen to illustrate our machinery 
on what may be one of the simplest examples of such an equation, the generalized 
parabolic Anderson equation, (gPAM), that was already handled in the 2-dimensional 


^We refer the reader to a more recent work by the authors [7] where such time-space paraproducts are 
more crucial to deal with higher dimensional space. More details are also provided there. Indeed, they 
consist of tensorial product between the semigroup and some standard 1-dimensional approximations 
of the unity to localize the time variable. 









41 


torus both by Hairer in [35] using his theory of regularity structures, and by Gubinelli, 
Imkeller and Perkowski in |32| , using their Fourier-based paracontrolled approach. This 
choice is motivated by the fact that only one (probabilistic) renormalization is needed to 
implement the paracontrolled machinery, while further renormalizations are needed in 
the stochastic quantization or KPZ equations. So the reader can see in the next section 
the machinery at work without being overwhelmed by side probabilistic matters. It 
makes sense now to make the following definition, in the present setting, where a 
stands for a real number in (0,1). 

► Definition. An enhanced, or rough, distribution C is a pair with C G ^ and 


The (generalized) parabolic Anderson Model in dimension 2 

This section is devoted to the study in our abstract setting of the (generalized) 
parabolic Anderson Model, in dimension 2. The analytical/geometrical setting is 
described in Section [2j The space (M, d, fj.) is a space of homogeneous type, equipped 
with a semigroup satisfying the regularity assumptions (lUEp and ( |Lipp . Let 

us insist here on the fact that even in this modest setting, the above semigroup approach 
offers some results that seem to be beyond the present scope of the theory of regularity 
structures, in so far as we are for instance allowed to work in various underlying spaces 
and even in the Euclidean space with operators L of the form div(AV), with A Holder 
continuous - see example 2 in section 12.11 The first two subsections are dedicated 
to proving some local and global in time well-posedness results, for the deterministic 
(gPAM) and (PAM) equations respectively. To turn that machinery into an efficient 
tool for investigating stochastic PDEs in which the singular term involves a Gaussian 
noise, we need to lift this noise into an enhanced/rough distribution; this step requires 
a probabilistic limit procedure generically called a renormalization step. It is performed 
in section 15.31 in the geometric framework of a potentially unbounded manifold, when 
working with a weighted noise. 

5.1 Local well-posedness result for generalized A big enough parameter b is fixed for 
PAM all the previous results to hold. 

► Theorem 5.1. Let a G (|)l) be given, and a' < a be close enough to a to have 
2q; -|- a' — 2 > 0; let also a finite positive time horizon T be given. Fix an initial data 
uo G and a nonlinearity E G (^^(M,]^). Given C, G set Z := TZ{0, and assume 

there is a sequence of smooth functions Cn converging to C in and a sequence 

of smooths functions Zn converging to Z in such that I[^^\Zn,Qn) converges in 

to some limit element Write Q = foi" tbe associated enhanced 

distribution. Then the generalized parabolic Anderson model equation (gPAM) 

dtu + Lu = (E(u), E'(u)tt') • C, tt(0) = Uo 

has a unique paracontrolled solution (u,u') G C" jq.^j(Z), with u' = F{u), provided T is 
small enough. 

Recall that the very notion of product in the right hand side of the (gPAM)equation 
depends on (^. Since, we have established in the previous sections the main analytic 
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estimates of paracontrolled calculus, we can prove this theorem l5.1l following the pattern 
of proof devised by Gubinelli, Imkeller and Perkowski in their seminal work [32], such 
as extended here to our more abstract setting. 


Proof — Given a singular spatial perturbation C € we know, by Proposition 

13.81 that Z := TZ{0 S Let 2l[o,r](Lii) be the set of controlled distributions 

{u,u') G snch that 

||(u,u')|L„' < K, u'{0) = F(uo) and (u,u')**(0) = uq. 

^a,[0,T] 

Denoting by v the solution of the equation 

dtv + Lv = (F(u), F'(u)u') • C, ^^(0) = uq, 
we define on 2l[o,r] (^) a map F setting 

^((u,u')) = (v,F(u)); 

it also takes its values in C" jq^j(Z), by Theorems 14.11 14.21 and 14.31 and has a 
C“^[o jn](.^)-norm bounded above by 


' "^[0,T] 

< |ho||c2c 


2) + 


{F{u),g 






+ F 2 ( ||F||c 3 ||u||g(<») (1 + ||C||c“-2) + I|F|Ic 3(1 + ||(■“J'“0||c“' )(l + ll■^llcTC“)(*) 


'[O.T] 


[0,T] 


with 

W := llCllc“-2 + l|F||c 3 |kllcTC“||Cllc “-2 + ||C^^^||c-rC2“-2' 

from the estimates provided in the above-mentioned theorems, together with Re- 

ol — ol' 

mark 13.101 describing how we can obtain the extra factor in terms of T 2 . The 
map is then a contraction from 2t[o,T](-^) to itself if K is chosen large enough 
and T small enough. Moreover, the Lipchitz estimates in Theorems 14.11 and 14.21 
give us the estimate 


d(3c,'[o u'), ^{v, v')^ 


< T- 



a,[0,T] 


((u,n'), {v,v')) 


for all {u,u'), {v,v') in with implicit constants depending only on K, Z 

and C- So $ happens to be a contraction of provided the time horizon T 

is small enough, from which it follows that has a unique fixed point in 2l[o,T](-fi^)- 
IFinally, it is easy to proceed as in Gubinelli, Imkeller, Perkowski’s work [32] and 
check that a solution of (gPAM) has to be in at least for a small enough 

T. 

> 


5.2 Global well-posedness result for linear PAM We focus in this subsection on the lin¬ 
ear (PAM) equation and prove a global 
in time well-posedness result in that setting. With that aim in mind, we define a 
weighted (in time) version of the previons functional spaces of paracontrolled distribu¬ 
tions. 
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Definition. Given A > 1 and a G (—2,2), set 

:= {/ G II/IIw.c- := sup e""* ||/(t)||^„ < ooj 

I t>o J 


and 

:= j/GLg'.C^'), 

We then define the space 


« := sup e 

s/t 

0<s,t<l 


-A. l|/W-/(s)|lc 


\t — s| 


< oo 


:= WAC“nwf L°°, 
;(«) 


and, given a reference distribution Z G VV’|“^ we define accordingly the space C^.^(Z) of 
pairs of distributions (/,/') in x such that 

V' 


(/,/')“ := / - nJ;^(Z) G Ct(WaC“+^) n 


Following the reasoning of Theorem 15.11 we prove in this section a global in time 
well-posedness result. One of the main ingredients used in the proof of Theorem 15.11 
was the Schauder estimates, through Proposition 13.81 or Corollary 13.91 We now give 
an extension of these estimates to the setting provided by the above exponentially 
weighted spaces. 


► Proposition 5.2. Consider an integer 6 > 0, /3 G (—2,0) and A > 1. For every e G (0,1) 
and V G W\C^ the function TZ{v)t := Pl;^gV{s) ds belongs to and satisfies 

the A-uniform bounds 

ll^(^)|UAC/3+2-2e ~ \HwxCP- 

and 

||^(''^)|| <3+^ 2^ < llullw c/3. 

Consequently, 

||7^(^;)|| 0+2 -2e < A~^ ||^|lw;,c/3- 

"w, 2 


Proof - We adapt the proof of Proposition 13.81 and add an extra new argument to 
consider the exponential weight in time. So consider another integer a > + 1, 

and a parameter r G (0,1]. Then 

Q(“)(7^(u)i) = fQ(^'>pEvis)ds. 

Jo 

Hence, 


g(“)p£u(s) 


< 

rs_/ 


< 


(—^ 

\T + t - S 

(—^ 

\T + t - S 


Q^it-sVis) 

OO 

(r + t-s)2||u(s)||^^. 


So by integrating, it comes 



s)2 dsj sup e ^^||'(^(.s) 11^,3 

/ sS[0,t] 

■S)2 ds) ||u||y^;^c/5■ 
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Let us just consider the integral term, temporarily denoted by I. If t < r, then 


I <t 2 f ds < r 2 L_^ < 

Jo 


< r2+l-^A-^ 


If t > r then 


I <t'^ e ds + J e 




^ ll — e 

< T2 -^- 


3 —Ar 


Jo 


^ ll — e £_i_i_f, _f 

< r2---h r2+^ ^A 

A 


e ^ dx 


< r2+1-^A-". 

So in both sitnations, we deduce that uniformly in A > 1 and t > 0, it comes 
,-A, qW(k(„),) < Ti+‘-'A-'||t,|U,c. 


and similarly 


-.—At 


Pi{nv)t) <A-^iiuiu,c/^. 


Consequently, we deduce that for every t > 0 

e"^1|^(^^)t|lc/9+2-2e ^ >^~"\HwxC/^^ 

which yields 

II'^^^)|IwaC/5+2-2s ~ ^W'^WwxCP- 

For the second statement, for s < t we have 

n{v)t - n{v)s = (p/!l - Id) {niv)s) + f P?X (x(r)) dr 

J s 

= j Q^^'^'R{v)s — + [ PfX{v{r))dr. 
Jo ^ Js 


So 


-.—At 


rQmv)s- <([ 

Jo ^ oo \Jo 


<1 I /2+^ — ]e ^‘||P(x)4c/3+2 


<{t- s)^+\X{t - s))"^||P(u)||^^^^+2 


{t - s)2+^-^X-^ 


WaC/5+2 
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and also (since /3 < 0) 




[ P^%{v{r))dr < j 

Js L°° Js 


< r Q^^\{r) — + 

't-r °o T 


P? (i;(r)) 


< 


< 


< 


At 


IvV’aC'5 


I dr 
r2-h 

T 


’Wllcn 

A't 

Js \Jt-r T 


P? (^^(0) 


+ l] dr 


dr 

dr 


IwaC/^ ^ e ((t - r)2 + l) dr 


< 




#+l-e 


where we used the fact that ^ + 1 G (0,1). So we conclude that 

||^(^^)|| P+i-e ^ ^"^II^IIwaC/9- 

w7 L°° 


l> 


► Theorem 5.3. Let a G (I , l) be given, and a' < a be close enough to a to have 
2 q; + a' — 2 > 0; fix also an initial data uq G together with some exponent A > 1. 
Given G set Z := TZ{C), and assume there is a sequence of smooth functions 

Cn converging to ( in and a sequence of smooths functions Zn converging to Z in 

such that n(^)(Z„,Cn) converges in WaC^" ^ to some limit element Write 
C = foi' the associated enhanced distribution. Then, for a large enough choice of 

parameter A, the linear (PAM) equation 

dtu + Lu = {u, u') ■ C, u(0) = no, 
has a unique global in time solution (n, n') G CP.^{Z), with n' = u. 


Proof - Since C £ C“ Proposition 15.21 implies that Z = TZ{C,) G Dehne 2tA(Ar) 

as the set of distributions (n,n') G C^.^^Z) with 

||(n,n')|Lc«' < K, n'(0) = no, (n,n')**(0) = no, 


and dehne $a((^,^‘ 0) ■= (^!^)) for (n,n') G 21a( iG), as the solution of v of the 
equation 

dtv + Lv = (n, u) ■ C, n(0) = uq- 

Then by the same considerations such as those detailled in the proof of Theorem l5.ll 
with some minor simplihcations getting around the non-necessary paralinearization 
step, and using Proposition 15.21 instead of Proposition 1,3.81 we get the following 
result. The pair ^\(u,u') belongs to C“.;^(Z) and satishes 


|$A(^^,^^')|| 


< 







with 

(*) := (llCllc“-2 + l|■^l|£AC“IICIIc“-2 + ||C^^^||c.j,c 2 “- 2 ) ’ 

uniformly in A > 1. A large enough choice of constants K,A ensures as a conse¬ 
quence that <1>A sends 21 a( iG) into itself. Moreover, we also have for (u,u') and 
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(v,v') in 'Qix{K) the Lipschitz estimate 


with implicit constants depending only on K, Z and C- So the result follows from 
the fact that <hA happens to be a contraction of ^x{K) if A is chosen large enough. 

[> 


5.3 Renormalization for a weighted noise We cannot expect to work in the Besov spaces 

used above when working in unbounded am- 
biant spaces and with a spatial white noise; so weights need to be introduced, with a 
choice to be made. We can either put the weight in the Holder spaces and still consider 
a uniform white noise, or we can put the weight on the noise and consider a weighted 
noise with values in unweighted Holder spaces. The first approach has been recently 
implemented by Hairer and Labbe in a forthcoming work on the linear (PAM) equation 
in see [38]. We choose to work with the second option here, partly motivated by ex¬ 
ploring this unexplored question, partly because it seems to us that spatial white noise 
in an unbounded space has more something of a mathematical abstraction than of a 
model for real-life phenomena. Refer to our follow up work [7| for the use of weighted 
Holder spaces in a paracontrolled setting. 

► Definition. Let a; be an L‘^{h) weight on M; the noise with weight u is the centered 
Gaussian process ^ indexed by such that for every continuous function / e 

we have 

(5.1) IE[^(/)2] = J f{x)uj{x)n{dx). 


Let us define the following notation. For t > 0, we denote by Qt the Gaussian kernel 


Qt{x,y) := 


1 


. d{^,vr 


V{x,x/t) 


it also depends on the positive constant c, although we do not mention it in the notation 
for convenience. Since we shall need to commute in some sense the Gaussian kernels 
with the weight lo, it seems natural to make the following assumption. We assume the 
existence of some implicit constants such that we have, for every t G (0,1] and every 
x,y G M, 


(5.2) uj{x)Gt{x,y) <uj{y)Gt{y,x), 

where we allow the implicit constants (in the exponential in Gt) to be different from 
the left to the right hand side, but but that they are independent of t, x, y. 

Recall the definition of Ahlfors regularity of a measure /x on a metric space (M, d ), 
given in section before the Sobolev embedding theorem 12.61 and quantified in equa¬ 
tion (|2.10l) . In that setting, it is relatively elementary to use the latter and prove by 
classical means that a weighted noise, as defined above, has a realization that takes 
almost-surely its values in some Holder space. 


► Proposition 5.4. Assume that (M, d,/x) is Ahlfors regular and let ^ be a noise on M, with 
weight Lo G nL°° satisfying the assumption (|5.2p . Then, for every cr < —there exists 
a version of still denoted by the same symbol, which takes almost surely its values in . 
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Proof — It suffices from general principles and lemma 12^ to check that the two expec¬ 
tations 


and 


® f / I® 

\jR‘i 


W :=IE( 


t 


are finite for every p > 2. We show how to deal with the second expectation, the 
first one being easier to treat with similar arguments. Starting from the fact that 
is, for every x G M, a Gaussian random variable with covariance the 
L^(u;^)-norm of K (a){x, ■), the equivalence of Gaussian moments with ()5.ip give 
the upper bound 


E 




(x) 


< E 




(x) 


< 

('N_/ 


J d<Qia){x,zfuj{z)p{dz) 


Using the Gaussian bounds for the kernel of with property (j5.2p and Ahlfors 
regularity, this implies for t G (0,1) 


E 


< w(x)2 t 


Hence, it follows that 


t P2< IIwIIe 


0 JM 


t 


1 P 

if cj < —the conclusion follows since w G n L°° C L^. 


U 


Let ^ be a weighted noise, with weight u, and define for every s > 0, a function 
5 s : M —)• M, by the formula 


so that we formally have 


gs{x) :=E[n(e-^^C,C)(x) 


roo p 

/ gs{x)ds = 1K Il{L~^^,^){x) 

Jo 


An explicit computation can be used in the case of the torus and the white noise to 
show that this integral diverges; see [35]. A similar computation can be done in our 
setting with the help of a highly non-trivial estimate on the kernel of the operators 
showing that the above integral also diverges at almost all points x of M. These 
facts justifies that we consider the modified integral (|5.3p below. Even though we shall 
only use here theorem 15.51 in a 2-dimensional setting, we prove it in the optimal range 
of homogeneous dimensions d G [2,4), for use in forthcoming works. Denote by S the 
function TZ{C) solution to the linear equation -|-A)H = ^, with null initial condition. 

Theorem 5.5 (Renormalization). Assume that (M, d,/r) is locally Ahlfors regular, with 
homogeneous dimension d = u £ [2,4). Consider ^ a weighted noise with weight a; G 
n L°° satisfying assumptions (15.2p . Set 


( 5 . 3 ) 


(H0C)(t) := f |n(^)(e ^'^^,C)(x) - 5^(x)| ds, 
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where we recall that gs{x) := E (e x) . Consider one of the following time 

functional space T = Ct, for an arbitrary finite time horizon T, or F = Wx, for some 
arbitrary A > 1. Then for every and p e (1, oo), we have 


E 




\p 

I J'C2“-2 


< OO. 


Moreover, by considering for e € (0,1), the regularized versions := e and := 
and 

poo 

:= I E ,f) ds, 


poo 

:= / E 
Jo L 

then for every p S [l,oo), we have 


lim E 
£—^0 




p 

jrc2a-2] 


= 0 . 


Note the following points before proceeding to the proof. 

• In particular, if the ambiant space M has hnite measure then the constant 
weight w = 1 satishes ()5.2H and belongs to r\L°°. So the previous results can 
be applied to white noise. 

• In Proposition 15.41 as well as in Theorem l5.5l we do not really need ui E 
it suffices that uj £ for a sufficiently large and finite exponent p. 

• Given a point o E M, any weight of the form uj{x) = (l + d{x,o)'j ^, satisfies 
the assumption (15.21) and belongs to L°° n , provided for M > |. 


Proof - By definition of white noise with weight w, we know that if T, T' are two 
self-adjoint operators then for every y,z £ M 


(5.4) 


E 


(TC)(y)(n)(z) 



Kt{u, y)KT'{u, z) uj{u)pL{du). 


Moreover if T and T' are self-adjoint operators, with a kernel pointwisely bounded 
by Gaussian kernels at scale t,t' S (0,1], meaning that we have 

\KT{x,y)\ <Gtix,y) and \KT'ix,y)\ < gf{x,y) 

for almost every x,y E M, then it follows from identity (|5.4p and Assumption (15.2p 
that we have 


IE[T^(y)T'^(z)] = J KT{u,y)KT'{u,z)u{u) p{du) < J Gt{u,y)Gt'{u,z)u{u) p{du) 
(5.5) <^{y) / Gt{y,u)Qt'{u,z) p{du) <uj{y)Qt+t’{y,z)-, 


we used Lemma lA.51 Fix now an integer 6 > 2 and for r E (0,1] and s > 0, 
define the quantity 
,"1 rl 

A{r, s) := 


0 Jo \r + tij \r + t 2 
tlt2 


Set 


{ti + s){t 2 + s) 

Qs :=n(')(e-^^e,e)-5s; 


, , x_4 dtl dt2 

{s + ti+t2) ^ir + ti+t2) 2———, 

t2 
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we claim that for every r G (0,1], s > 0 and every x £ M then we have 


(5.6) 


E 




< A{r,s)u}{x f. 


Step 1 — Proof of (|5.6I) . The resonant, or diagonal, part of the paraproduct 
is given by five terms, of the form 

nHf,g) = £vt{itL)Vlf-Qtg)j or nHf,g) = £vt{Qtf-mvh)j, 
or 

9) = £ VtT (VtVlf , VtVh) j 

where 

• Vt,Vl and Vf are operators of the form p{tL)e~^^ with p a polynomial func¬ 
tion; 

• Qt is of the form {tLY~^p{tL)e~^^ with a polynomial function p. 

So both of these operators have a kernel with Gaussian pointwise estimates and we 
only need to deal with these three generic quantities. 

Let us focus on a term of the first form and study 


e] :=n^{e - gl{x) with gl{x) :=E n^{e 


Due to the covariance rule of Gaussian variables, we have for T, U, T', U' self-adjoint 
operators (using (I5.4p i and every y,z £ M 


E 


Hence E 

1 rl 


'my)Uay)r^{z)U'az)-E[my)Uay)MT'az)U'az)] 

= E[TC{y)raz)]E[UC{y)U’az)] + E[rC(y)C/'e(^)]E[[/e(2/)n(^)]. 

|Qr^^0^(x)|^ is equal to 





0 JO 


where we use the notation /(•,•) for a function of two variables, with (fg)(»^») 
standing for the map {y,z) f{y,z)g{y,z). Moreover, to shorten notations, we 

shall use below the notation dm for the measure p{dy)p{dz)^^. By applying 
(ESD, it follows that 

| 2 ' 


E 


\Qi^^elix)f 


<Jl + J2 


with 


/ i 1 ^2 

w(j/)w(z) j-^^j-^^gt^+t^+s{y,z)Qt^+tAy^z)dm 


and 


/ i 1 ^2 

j-^j-^gt,+t2+s{y,z)GH+t2+s{y,z)dm. 
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Let us first explain how we can estimate the kernel of Using the notation 


Pfj = p{tiL)e for some polynomial function p, it comes 
Qr^'Pti = ((^ + e-^^p{tiL)t 


—t\L 


r + Ui 


(b) 


has a kernel with Gaussian bounds at the scale 

(b). 


so since r + ifi ~ r + U, Q 
r + \ti and p(fiL)e“ 2 *U scale fi, it follows by Lemma rA.5l that Qr^^Vti has 


r+ti 


a kernel pointwise bounded by Gr+ti with an extra factor 
to estimate the first term Ji. We have the upper bound for Ji 

^2 b 


Coming back 


(r + ti)(r + t2) 

<bo{xf j 


Gr+u (x, y) Gr+t^ (x, z) uj{y) uj{z) —— —— Gn+t^+siy, z) Gt^+t-, (y, z) dm 

+ S 12 + s 


(r + ti)(r + t2) 




^2 


■ Gr+ti {x, y) Gr+t2 {x, z) Gti+t2+a{y, z) Gti+t2 {y, z) dm, 


ti S t2 s 

where we used Assumption (2121 Due to Lemma [A.51 with Ahlfors regularity (|2.1U|) . 
we have 

jj Sr+ti {x, y)Gr+t 2 {x, z) ati+t 2 +s(y, z) Gti+t 2 {y, z) p{dy)p{dz) 

{ti +t 2 + s)~i JJ Gr+ti {x, y) Gr+t 2 (x, z) Gti+t 2 {y, z) p{dy)p{dz) 


< 


Hence, 


Ji ^^{xf 


1 rl 


^ (^1 + t 2 + s) 2 (r + ti + t 2 ) 2 • 
tl t2 


0 JO 


{r + ti){r + t 2 ) 


{s + ti+t2) 2(r + ti+f2) 2— — 
U + s ^2 + S ^ 1^2 


<u]{xfA{r,s). 

The second term J 2 can be similarly bounded, which concludes the proof of (I5.6D 
for 0^. The corresponding term 0^ with Ti? can be estimated in the same way. So 
it remains us now to focus on the last and third term with 

n\f,g)= fvttT{vlf,Vh) Y 

Jo b 

and 0^ := Tl^[e~^^^,^){x) — gs{x). Following the exact same reasoning we have 


E 


.)|") <^1 + ^:2 


with Ki equal to 


and K 2 equal to 




tit2 


{tl + s){t2 + s) 


tit2 


GH+t2+s{y, z) GH+t2{y, z) dp 


Gti+t2+s{y, z) Gt^+t2+s{y, z) dv. 


{tl + s){t2 + s) 

Following the same computations, gives us that both Ki and K 2 are bounded as 
follows 


All + A :2 < u}{x)‘^A{r,s), 
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which concludes the proof of (15.6p . 


Step 2 — Conclusion. We refer the reader to Lemma lA.61 for a precise control of 
quantity A. Combining (I5.6p with Lemma lA.61 gives 
(5.7) 


E 


<E <a;(x) 


s + r 



1 + log ( 



We then consider 

(HO?)(t):=y - gs{x)^ds = J Qsds. 

We refer the reader to Definition 12.41 for the definition of Besov spaces. For all 
0 < s < t, it comes (where we forget the low-frequency part in the Besov norm for 
simplicity, since it is easier than the high-frequency part which we deal with) 


E 


\Eom-^oas)\ 


2p 


B: 


•2a —2 

2p,2p 


= ; r 


-p(2a-2) 


E 


IM 


QW[EOe(t)-HOC(s)] g{dx) 


2p 


dr 


< 


j\-pi2a-2) J J^^*]E(^|QW0(r)(x)|) dr 


2p 



< 



^-p{2a-2) 




1 

2 



dr 

r 


where we have used Gaussian hypercontractivity, that is the basic fact that the 
LP norms are all equivalent in a fixed Gaussian chaos. So it comes, by Minkowski 
inequality. 


E 


|H0e(t)-H0e(s)||J.-2 

^2p,2p 


|2p 

l2p 


t / 


„-2p(a-l) 


r -|- r 


j_ <. 2p 

+ dr] , 


We have 

pi 


-2p{<y—l) 


/O 


< .j.-2p(a-l)-dp J 


P dr 
r 


r\P^ dr 
r 


< -2p(a-l)-dp 


where we chose an arbitrary small parameter e with 

—2(n — 1) — — -|-e<0< — 2(q; — 1) + 1 — — 


which is equivalent to 


, d e 3 

'-4+2<“<2- 


d 

4 


and such e > 0 exists. Observe that the latter condition is satisfied since d £ [2,4), 
sol — |<a<2 — |<| — |. Then because a < 2 — so —{a — 1) — | > —1, we 
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have in the end 


E 




2p 


< 


|2p| 

\2p\ 


|W||ot,|s — t\ 


—2p{a—2)—dp 


We can then use Kolmogorov’s continuity criterion to deduce that for every T < oo 
and A > 1, we have 


E 




2p 

CtB. 


2a-2 

2p,2p 


+ E 


l-oe| 


2p 

SxB. 


2a-2 

2p,2p 


< oo. 


As we are in a position to use Besov embedding, as we assume Ahlfors regularity 
()2.10p . we know that 


d2q-2 n2a-2 ^ 2p _ ^2a-2-^ 

-°2p,2p ^ '°2p,oo ^ -Doo,oo — O 


So for every p > 1 and every aG (l — ^,2 — 


E 




2p 

CtC^ 




+ E 


l-oe| 


2p 


■2p 


< oo, 


which allows us to conclude. 

Step 3 — Second part of the statement. The very same reasonings prove the 
second part of the statement about the approximation results. We only sketch the 
arguments. First, observe that estimating the main quantity 

is almost equivalent to the study of 

Then by the bilinear structure, this is equal to 
(5.8) (E-E")0e + H"0(C-f)- 

We may then repeat the two previous steps to estimate the two terms, indepen¬ 
dently from one another. Note that E — S'" is the resolution of the heat equation 
associated with ^ so both of the two terms make appear the quantity ^ 

So to apply the two previous steps to any of the two terms in ()5.8p is very similar 
and we only have to include some additional factors coming from 

da 


(5.9) 


e-r = (l-e-"^)C= [ 

Jo 


a 


A careful examination shows that if we replace in the previous reasoning (steps 1 
and 2) the terms ^ by then all the estimates are identical and we get an 

extra factor of the forms 

sa ati at2 

(s-|-cr)2’ {a + tiY ° {<y + t 2 Y 

So in these three situations (by replacing ^ with the same estimates hold 

with a quantity A(j(r, s) satisfying, for rj > 0 as small as we want the inequality 


(5.10) 


Aa{r,s) 


s + r 


(rs) 2 I 1 + log 


s + r 


a 


mmir, s 


instead of estimate (jA.9p . Such an estimate can be plugged into p5.9p and provides 
an integrable quantity on an interval r G (0,e); this gives similar estimates for 
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as those obtained for H<C>^, with an extra factor in terms 
of and r], as small as we want, is also appearing in the conditions on a,p ,.... We 
then let the reader check that since all the conditions on the exponents are “open 
conditions”, then the previous reasoning can be reproduced, up to a small loss of 
regularity. That means that we can incorporate such a sufficiently small parameter 
rj > 0, and with 

we can get 


E 




2p 

B: 


•la — 2 

2p,2p 


< 



0 Js 




2p 


a J 

2pri 


which allows us to conclude as previously using Besov embedding. 


O 


Proof of Theorem 11.21 - Fix the weighted white noise ^ and its regularized version 
:= As in Theorem 15.51 or Theorem 11.21 set 


:= 


E [n(e-"^f)] ds = E[n(L-if,f) 


In order to make appear this term in the equation, we can introduce a suitable 
correction term in the regularized problems, which leads us to study the following 
renormalized PDF 


dtu^ + Lu^ = - C^F'{u^)F{u^); 


this is elementary, and detailed in [32] . We can then follow the exact same approach 
as for Theorem 15.11 or Theorem 15.31 for the global estimates with the spaces £\, 
adapted to this modified PDF. So we only detail the required modification. We 
cannot use which does not converge, and use instead — 

which converges in 

Using the arguments of Theorem 14.11 with 7 = a — 2, one defines a new “product” 
for (u,u') e C“(A") 


((F(u),F'(u)u') - f) -CVF'Cu) G 

(whose the norms can be estimated by quantities independent with respecto e), for 
which we have the uniform estimate 


{F{u),F'{u)u') ■ e - n^(„)(r) - C^u'F'iu) 




< 


{F{u),F\u)u') (\\e\\c^-^ + \\X^c.c^\\e\\c^-^+ V 

We conclude as in the proofs of Theorems 15.II and 15.31 by a fixed point argument. 

U 
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A _ 

Heat kernel and technical estimates 

We gather in this Appendix a number of propositions whose proofs were not given 
in the course of the paper, so as to keep focused on the most essential aspects of our 
work. These proofs are given here. 

We start by proving the following pointwise and L^-estimate for the gradient of the 
heat semigroup. 


Proposition A.l. Assume that (M, d,//) is a doubling space equipped with a semigroup 
satisfying (lUEp and ( Lip ). Then for every t > 0, xq £ M and every function / e we 
have 
(A.l) 


{Vir){e *^/)(xo) 



1 


V[xQ,Vi)V[y,y/t) 


exp 


—c ■ 


d{xo,yf 


\f{y)\n{dy). 


Let us first introduce the following notation: for a function / £ and a ball 
B C M, we write Oscs(/) for the oscillation of / on i? dehned by 


Oscs(/) := f-J^fdn 


dy, 


where f dy stands for the average of / on the ball B. 

Proof — Fix the function f £ and consider g = e“*^/. By L^-Caccioppoli inequality 
(see Lemma below), we have for every xq and r > 0 that 


/ r(5)2(i/r) <-OscB{xo,2r){9)+(-f \L9\ 

J B{xQ,r) j ^ \J B(xn,2r) 


B(xo,2r) 


dy 


So if xq is a Lebesgue point of r(( 7 )^ and \Lg\‘^ (which is the case for almost every 
point Xq £ M) then taking the limit for r —)• 0 yields 


(A.2) 


r(^)(xo) ^ lin^nf ^ 2r ){9) + \L9{xo)\. 


Since {tL)e has a kernel satisfying the Gaussian upper estimates (jUEp (by an- 
alyticity), we deduce that 


(A.3) \tLg{xo)\ < f 

^JV{xo,Vt 


Vi)V{y,Vi) 


: exp —c 


dixo,yf 


\f{y)\y{dy). 


Using Lipschitz regularity ( |Lip[ ) for the heat kernel and doubling property, it comes 
for x,z £ B{xo, 2r) with r < ^/i 


\g{x)-g{z)\= e *^/(x) - e ‘^/(z) 


< 

CN-/ 


d{x, z) 


^ d dM ^v{xo,Vi)V{y,Vi) 


exp ( ) \f{y)\y{dy). 
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Hence, uniformly with respect to r G (0, y/i) we obtain 
(A.4) 

^OscB(,„,2r)(5) ^ 4 / / . .=^ exp ^ \f{y)\^{dy). 

^ VtJM^V{xo,Vt)V{y,Vi) V i J 

By combining these last inequalities (lA.lih and (IA.4p into (|A.2|) . one concludes to 

mi- 

t> 

► Lemma A.2 (Cacciopoli inequality). For every ball B of radius r > 0 and every function 
/ G 'D 2 {L) we have 

j. 1 

(£ Tiff d^,y Osc 2 b(/) + r |L/|2 ' . 

Before to check this inequality, let us first recall some consequences of the Gaussian 
upper estimates (lUEp . Under (jUEp . we know that a scale-invariant local Sobolev 
inequality holds, more precisely 

ll/ll' <l^l"''(ll/ll2+^'^(/,/)), 

for every ball B of radius r > 0, every / G 'T) 2 iT) supported in B and for some q > 2. 
This inequality was introduced in [56] and was shown, under (IVDh . to be equivalent 
to (|UEp in the Riemannian setting. The equivalence was stated in our more general 
setting in |SS]. See also [TB] for many reformulations of local Sobolev inequalities, an 
alternative proof of the equivalence with (lUEp . and more references. 

Such a local Sobolev inequality also implies a following relative Eaber-Krahn inequal¬ 
ity (see for instance [421 Theorem 2.5], as well as [161 Section 3.3]): for every ball B 
with a small enough radius r > 0 , every function / G 'T>2{T) supported in B then 

(A.5) ll/lb <r||r(/)||2. 

Proof of Lemma IA.2I - We refer to m Lemma A.l] for such a result for harmonic 
function: if tt G 1*2 (A) is harmonic on 2B (which means Liu) = 0 on 2B) then 

(A.6) ^jTiu)‘^df?j <^Osc2b(u). 

Now consider / G P 2 (E). By [111 Lemma 4.6], it is known that there exists 
u G T) 2 iL) harmonic on 2B such that f — u & L^iV) is supported on the ball 2B. 
By the support property, it follows 

||r(/ - u)||2 = -u)dfJ, = J if - u)Lif) dfj,. 

So using Eaber-Krahn inequality (IA.5I) we obtain 
(A.7) 11/- “lb ^ ^||r(/- «)IL2(2 b) 

and so 

||r(/ - u)\\l < 11 / - «||2||a(/)||^3(2b) ^ H|r(/ - “)|| 2 IIa(/)IIlr 2 b), 

which yields 

(A.8) l|r(/-f)lL<-'||n/)lli.(2B)' 
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Then we split 

||r(/)|i, < ||r(/- .<)||, + lirwii, < + llrWlL 

and then use (lA.Bp to get 

< r (j-^^\L{f)\^ df?j + ^ 0sc2b(/) + ^ \f - df?j 
1 

^ l^(/)I^C^/^) +^0 sC2b(/), 

where we used again (jA.7p and (|A.8p at the last step. 


O 


We also give a proof of the following basic important fact about the Holder spaces 

Proposition A. 3. For cr < 2, the Flolder spaces do not depend on the parameter o 
used to define them, and the two norms on corresponding to two different parameters 
a, a' , are equivalent. 

Proof - Given two positive integers a and a', consider the two spaces and and 
their corresponding norms. Fix t S (0,1]. If a' > a, then writing 

2 2 2 

and using the fact that the operators Q^t are uniformly bounded on L°°, we 

2 

get 


\r’a • Wna 

|G , rv_/ II IIL-ct 


If now a' < write 


{a') _ 1 


QT’ = 


la—a' Jo 


g(a')g(a-a') ^ ^ . 


For s < t, we have 


K)^(a-a') ^ /'sy “ g(a) f_J_ 


Qrvi 


t-\-s 


t S 


so that for f € 


gK)g(a-a')^ 


< 

r\j 


< 


Qltf 

t^WfWc,. 


For t < s, we have 






t s 
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SO that 


and similarly 




< 

oo V 5 


< ( - 
5 


Qtlf 


S2 


gia') p{a-a') J, 


< 


Then by integrating (and since a', a — a' > 1 > f) we have 




fi 




t\ ^ ds 

- S2 - 

t \S) s 




^ iiJiicg 

which concludes the proof that 


h»cr 

^ / rs_/ 

n' 


\> 


The following lemma provides a useful way of proving that a distribution is Holder; 
it was used in sections 12.31 and 15.31 to investigate the almost sure regularity properties 
of white noise and the renormalized paraproduct dealt with in theorem [53J Recall that 
Besov spaces were defined in Definition 12.41 

Lemma A. 4. Assume that the metric measure space (M, d,/r) is Ahifors regular (see 
()2.inp h with exponent v. Then, given —oo < cr < 2, and 1 < p < oo, we have the 
continuous embeddings 


B, 


p,p 


^P,00 ^ r>00,CXD — O 


Proof — The first embedding is a direct application of the following fact. For s G (0,1) 
and an integer a >2 then 

Qi“V = - 

S 

Since the semigroup is uniformly bounded on U’, we get 




< 

rsj 




dt 

p t 


and by Holder inequality 




< 

rs_/ 




P dt 
p t 


< 


S2 


■ 

P-P 


The second embedding comes from the following elementary fact. For t G (0,1), 
let T a linear operator with a kernel, pointwisely bounded by a Gaussian kernel Qt 
at scale t, then with Ahifors regularity (I2.10jl . we have 


rii 


LP->L° 


<t 2 p. 


So for s G (O, ^), applying to T = we obtain since 




< Q 2p 




LP 


_I 

< S 2p 2 

r\j 


^p,oo ‘ 
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which proves the embedding oo 


B 



C0,00‘ 


!> 


The next three lemmas were used in the proof of the renormalization theorem 15.5 


Lemma A. 5. For every t > 0, set Qt the Gaussian kernel at scale t 


Qt{x,y) := 


1 


V{x,Vt) 


d(^,vr 

e ‘ , 


where we forget the dependence with respect to the constant c in the notation. Then for 
s > t > 0 and every x,z £ M, we have 


Qs{x,y)Qt{y,z) y{dy) <Qs(x,z) 


< 


1 


V(x,^} + V(z,./s}' 


Proof — By considering that Qt{-,z) belongs to uniformly in s, we directly obtain 
that 

1 


Qs{x,y)Qt{y,z) y{dy) < sup Gs{x,y) 


< 


y&M 


V{x,y/s)' 


Moreover, 


, d{x,yf\ ( d{y,zf\ ^ ( d{x,yf\ ( d{y,z^‘^ 

exp ( --- ) • exp (- - - I < exp (- I • exp ' 

d{x, zY 
2s 

So in the product Qs{x,y)Qt{y-, z), we may factorize an exponential decay and so 
for some implicit constants, we have 


< exp — 


J Gs{x,y)Gt{y,z) ix[dy) <Qs{x,z). 


o 


Lemma A. 6. For r G (0,1], s > 0 and d> 2, let us consider the quantity 
»1 />! 

A{r^ s) := 



/o do \x + tij \r + t2 
tlt2 


V (^1 + s){t 2 + s) 


(s + + t 2 ') 2 (r + + t 2 ) 2 


where 6 > 2 + | is an integer. Then we have 


(A.9) 




dti dt2 
ti t2 ’ 


Proof - The two variables ti , t 2 play a symmetric role so we may restrict our attention 
to the double integral under the condition t 2 < ti- The part Ai of the double 
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integral where t 2 < ti < r gives 
/ ht2 


Ai(r,s) = 



10 


0 V (^1 + ®)(^2 + s) 


2 _d d dt2 dti 


< 

rsj 


< 



ti 


ti + s 


(tl + s) 2r 2 


tl 


ti + s 


l + log(^)j 


s + tis \ dti 


tl 


I' ( (t.+VO ('+“■'5 (^)) 


< 


2 |s +log ('S+ 


< (sr) 2 


s + r 

where we used the basic inequality 


I 


dt2 


< 


0 y/t2{t2 + s) \S + tl 


tl 


1 + 


>-(^)) 


which can be easily checked by splitting into the two cases ti < s and s < ti. The 
second part A 2 of the double integral where t 2 <r <ti \s controlled as follows 


M{r, 



t2 h 

/ -4 dti 

ih + s) — 

tl 


< 


r + s 


(r(r + s)) 2 + log 


where we used Lemma [A. 71 The third and last part A 3 of the double integral where 
r < t 2 < tl satisfies 


^3{r,s) 



tlt2 


{tl + s){t 2 + s) 
r 


-f,, , ._ddt2dti 

h^{ti + s) 2 __ 
t2 tl 


r + s 


, .-d -4 dti 

{ti+s) . 

tl 


We then use again Lemma lA.71 to obtain 


^3 (as) 


< 


r + s 


(rs)- 


Inequality (IA.9P comes by combining the above three estimates. 


O 


Lemma A.7. For every r, t G (0,1) and any 0 < p < e 



t -\- S 


^ sdt ^ 


t 


s + r 


Proof — Indeed if r > s then for every t G (r, 1) we have t ~ t + s and so 

"1 / + \p 



t 


t s 




< f\-^— 

t ^ Jr t 
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Now if s > r, we split the integral in two terms and we have 



where we used e > Ap. 


I> 


B 


Extension of the theory 


Consider as above a doubling metric measure space (M, d, p) equipped with a heat 
semigroup satisfying the upper gaussian estimates (|UEp . We aim in this appendix at 
explaining how one can get the same conclusions as in the above main body of work 


(a) by weakening the Lipschitz regularity assumption on the heat kernel (Lip), 
assuming only some integrated estimates of the gradient of the heat kernel; 


(b) by developing the theory of paracontrolled calculus in Sobolev spaces rather 
than in Holder spaces. By Sobolev embedding, Sobolev spaces are included in 
some Holder spaces, so it will be interesting to understand if starting from an 
initial data belonging to some Sobolev space, the solution of renormalized singu¬ 
lar PDEs will lives in this same scale of Sobolev spaces. From a technical point 
of view, it is a bit more difficult since Sobolev spaces involve simultaneously 
all the frequencies, whereas for Holder spaces we can work at a fixed frequency 
scale. 


We give in section IB. II the regularity assumptions on the heat kernel under which we 
shall work here, and reformulate and extend in section [B.2l the main continuity estimates 
on the operators and E needed to extend the paraproduct machinery to 

the present setting. The latter, together with some crucial commutator estimates in 
Holder and Sobolev spaces, is investigated in section IB. 31 The last and short section 
IB.41 describes how these results can be used to extend the results of section [S] to our 
optimal regularity setting. 

This appendix was written jointly with Dorothee Prey. 


B.l Regularity assumptions Rather than assuming the Lipschitz property (Lip) used 

above we shall assume here that the gradient /carre du 
champ operator T satishes some estimates and the L'^-de Giorgi property recalled 
below in sections IB. 1.1 1 and IB. 1.21 We shall also assume that it satisfies a scale-invariant 
Poincare inequality recalled in section IB. 1.31 

B.1.1 L'^-estimates of the gradient of the Given qq > 2, the uniform -boundedness 
semigroup of the gradient (or “carre du champ”) of the 

semigroup was introduced in [2] 

(Gpo) sup||Vtre"*-^||go^go < -hoo. 

t>o 

By definition of the carre du champ operator, (G 2 ) holds trivially. It is known in 
that case that this global L^-inequality can be improved into localized estimates, via 
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L^-Davies-Gaffney estimates. For every subset E,F C M and every t > 0, we have 


^-tL\ 




+ v^llre 


-tL\ 




< e"'^ 

r\j 


cP{E,F) 
t 


Assuming the volume doubling condition (VD) and the Gaussian upper bound (UE) 
for the heat kernel, one can interpolate the estimate (G^p) with the above L^-Davies- 
Gaffney estimates and deduce that (Gg) holds for every q G [2, go]- More precisely, for 
every subset E, F C M and every t > 0, the ineqnality 


holds for some positive constant Cg, only depending on g e [2,(?o)- Following [21 Propo¬ 
sition 1.10], the latter estimate can be reformulated in terms of integral estimates of 
the gradient of the heat kernel. Denoting by pt the kernel of we have 


Vi\\r^pt{.,y)\\^ < 


V (y, Vt) 


for /r-almost all ?/ G M and all positive times. By interpolation with the L^-Gaffney 
estimates, there exists a positive constant c such that 


(B.l) 



„ d(-,i/)2 


^xPt{;y) 


< 

rN_/ 


V (y, Vi) 


q L j 

holds for //-almost all y G M and all positive times. 

We refer the reader to |2] for more details about Property (Gp) and the link with 
the boundedness of the Riesz transform; see also m and references therein for more 
details. 


B.l.2 L^-de Giorgi Property The so-called ”de Giorgi property”, or ’’Dirichlet property”, 

on the growth of the Dirichlet integral for harmonic functions 
was introduced by De Giorgi in [22], for second order divergence form differential oper¬ 
ators on M”’, with real coefficients. In de Giorgi’s work, this property prescribes a(n at 
most) linear growth rate for the L^-average of gradients of harmonic functions. This 
property was subsequently used in many works and in various situations in order to 
prove Holder regularity for solutions of inhomogeneous elliptic equations and systems. 
An L'^-version was recently introduced in nn, and we refer the reader to that work for 
more details about it. 


Definition (L'^-de Giorgi property). Given q G [l,-|-oo) and 9 G (0,1), we say that the 
operator T satisfies the inequality (TDG^)^ if it satisfies the following estimate. For every 
positive r < R, every pair of concentric balls with radii r and R, respectively, and 

for every function f GV, one has 

iDGg,e) ~ ( 7 ) {(£ |r/rd/^)'+i?||L/|Uoo(s^)|. 

We sometimes omit the parameter 6, and write {DGg) if (IDGg ^l) is satisfied for some 
0G(O,1). 


As we always have 

for every f G R and 0 < r < ii, if the space is doubling, with dimension ly, the 
inequality {DGgfi) holds for every q> v, with 9 = ^ < 1. 
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B.1.3 Poincare inequality 


(^2) 


Last, we shall assume that the carre du champ operator F sat- 
ishes the following scale-invariant Poincare inequality 


L 


f-i 




< 


'B 




for every / G 2^2(L) and every ball B of radius r. We refer the reader to |llj for 
a precise study of the connection between Poincare inequality, L^-gradient estimates 
and de Giorgi property. Let us just point out that if P satisfies the above Poincare 
inequality and the gradient estimate (Gpu), then there exists a parameter 6 G (0,1) 
such that the inequality (DGp^e) holds for every p G [2,po). Note also that in the first 
and main part of this work, we assumed an upper Gaussian pointwise estimates for the 
gradient of the heat kernel equivalent to (Goo)- This assumption yields the Poincare 
inequality {P 2 ), the integrated gradient estimate (Ggg) for every qq G [2, 00 ], and also 
de Giorgi property (DGq^e) for every 9 G (0, 1) and every q G [2, 00 ). 

Our aim in this appendix will thus be to weaken the (Goo) assumption made above 
into a combination of (Ggp) and {DGq^g), for some exponent qQ,q and 6. 


In the first and main part of this work, the paracontrolled calculus and its applica¬ 
tion to the 2-dimensional parabolic Anderson model equation was studied under the 
assumptions (lUEp and (Goo) that the heat kernel and its gradient satisfy pointwise 
Gaussian upper bounds. The aim of this section is to weaken the latter condition. 
Here are examples where the operator P satisfies only the properties {P 2 ), (Gq) and 
(DGq) for some q > 2, and where (Goo) does not hold. 

(a) Conical manifolds. Gonsider a compact Riemannian manifold N of dimension 
n — 1 > 1, and define M := (0, 00 ) x N as the conical manifold whose basis is N. 
It is known that M is a doubling manifold of dimension n which satisfies (jUEp . 
Moreover, as shown by Li in [37], the operator P satisfies (Gg) if and only if 


q < q{N) := 


n 


n 


where Ai is the first non-vanishing eigenvalue of the Laplace operator on N. As 
an example, if we consider N = rS^ the circle of radius r > 1, then 

2r 


q{N) = 


r — 1 


(b) 


So theorem IB. 181 below allows us to solve the PAM equation on M for r suffi¬ 
ciently close to 1. 

Elliptic perturbation of the Laplacien. On the Euclidean space or 
any non-compact doubling Riemannian manifold satisfying Poincare inequality 
(P 2 ) and the Gaussian bound (|UEp . we may consider a second order divergence 
form operator L = —div(AV) given by a map A taking values in real symmetric 
matrices and satisfying the usual ellipticity condition. Then if A is Holder 
continuous, it is known that —L generates a self-adjoint semigroup with (|UEp 
and Gaussian pointwise bounds for the gradient of the semigroup (Goo); see 
H]. In such a case we may apply the results proved in the first part. Eollowing 
Auscher’s work [T], we know that the combination of property (lUEp with Holder 
regularity of the heat kernel is stable under perturbation. So fix Aq a 
Holder continuous map with values in real symmetric matrices and satisfying 
usual ellipticity condition. Then for every Q > 2, and any positive 0, there 
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exists a positive constant e such that for any map A on the state space, with 
values in the space of real symmetric matrices, and such that ||^ — Ao||oo < 
the operator L = —div(AV) satisfies (Gg^) for some go > Q, and has de Giorgi 
property (DGg^e), ior 9 = d/q < Q and some Q < q < go- In such a situation, 
we may apply Theorem IB.181 and deduce that we can solve the PAM equation 
in such a 2-dimensional context. 

(c) Lipschitz domain with Neumann boundary conditions. Similarly, con¬ 
sider an open and bounded subset C and consider for L the self-adjoint 
Laplace operator associated with Neumann boundary conditions. Then by a 
change of variable, this situation is very similar to the previous one: if the 
boundary is sufficiently close (in a Lipschitz sense) to a smooth set (at least of 
regularity C^), then we can solve the PAM equation. 


B.2 Functional calculus and gradient estimates in We start this section by quantifying 
Holder and Sobolev spaces the localization properties of the op¬ 

erators and , and their gra¬ 
dients, in IjP spaces, before turning to the gradient estimates of the heat semigroup in 
the intrinsic Holder and Sobolev spaces in section [B.2. 21 


B.2.1 Localization properties of the ap- As we know, for every integer a > 0, the operators 

proximation operators and and have a kernel satisfying Gaussian 

estimates (|UED . The above regularity assump¬ 
tions (Gq), (DGq), {P2) on the gradient operator actually imply much more. 


► Lemma B.l. Let pt stands for the kernel of or P^°‘\q[°‘'^ for any integer a > 1. 
Under {DGq^e) and (Ggp) with Poincare inequality (P 2 ) foi' some 2 < q < go, we have the 
following Holder regularity estimate for the heat kernel. For every r/ S (0,1 — 0], f > 0 and 
almost every x,y,z ^ M 


\ptix,z) - ptiy,z)\ 


d{x,y) 

Vi 


Viz, Vi) 




d{x,z)'^ 


We only sketch the proof and refer the reader to mi for details. 

Proof — We follow the argument of Morrey’s inequality, which relies oscillation es¬ 
timates to some gradient bounds. Let x,y £ M be Lebesgue points for / = 
Pti-,z) with d{x,y) < Vi, otherwise there is nothing to be done. Let i?j(x) = 
H(x, 2“*(i(x, y)), for i £ N. Note that for all f £ N, i?j(x) C Bo(x). By Poincare’s 
inequality, this yields 


fix) - / fdp 

J Bo{x) 


^2 ^dix,y) (-f 

i>0 \JBi{x) 


\rf\Up 


By considering a ball of radius Vi containing both x, y, iDGq^g) yields 


V 

jBiix) 


\rf\Up 


< 


Vi 


2-^dix,y) 
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Since f = pt{-,z), by (Gg^) and (jUEp we know that 

1 

SO we can conclude the proof by summing over i, since 0 G ( 0 , 1 ). 


O 


Under the sole assumption (UE) that the kernels of the operators and have 
Gaussian upper bounds, these operators are bounded in every L'p space for p G [1, oo], 
uniformly with respect to t G (0,1]. Moreover, for every pi,p 2 £ [IjCc] and t > 0, 
they satisfy the following off-diagonal estimates at scale y/t, which quantify 

the localization properties of these operators. Eor every ball Pi,P 2 of radius y/l, and 
for every function / G we have 

(£ (£ ifrdpY' . 

One can refine this estimate by using off-diagonal estimates, such as done in m Lemma 
2.5, Lemma 2.6]. 


► Proposition B.2. Assume (Gg^) for some qo > 2. 

(i) For every non-negative integer a and every p G [2,go), the operators rp/“^ and 
satisfy L^-L^ the following off-diagonal estimates at the scale y/i. For every 
ball Bi,B 2 of radius y/i and every function / G L^(Pi), we have 


£ \yftrQ[^'>f\Pdpy + (^£ \yftr4^^f\Pdp 


< e-^ 


d(fll,g2) 


/ I/P 

Jb^ 


dp 


It follows in that we have 


sup 

t>o 


(Vtr)(p/“).) 


+ 


p^p 




p^p 


< oo 


for every p G [2, go)- 

(ii) For every positive real number a and every positive t, the operator is an integral 
operator with kernel that satisfies the inequality 


(B.2) 


kt\x,y) 


< 

rsj 


V{x,y/t)^V{y,y/t) 


1-X 


1 + 


d^{x,y) 


for all A G [0,1] and /i-almost all x,y £ M. As a consequence the operator 
satisfies the following LP^-LP'^ off-diagonal bounds of order a at scale y/i, for 
every pi,p 2 G [1,-Foo]. Given any balls Pi,P 2 of radius y/i, and any function 
/ G P^i(Pi), we have 


/ iQ!“Vr 

JB2 


dp 


^ I 1 + 


d(Pl,P2)' 



Besides these localization property in the physical space, the approximation opera¬ 
tors satisfy some orthogonality properties, which will be of crucial use in proving 
the continuity properties of the paraproduct and resonant operators below, and which 
can be viewed as an analog of the Littlewood-Paley theory, as made clear in [121 Propo¬ 
sition 2.13 and Lemma 2.15]. 
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► Lemma B.3. Let o be a positive real number. Set 

Qt := = 2tQf, 

2 

so that = Q^. Assume the Gaussian upper bound (lUEl) holds. Let also F : (0, +oo) x 
M —M be a measurable function and write Ft{x) for F{t, x). Then for every p G (1, +oo), 
one has 



whenever the right hand side has a meaning and is finite. If F = / does not depend on t, 
we have the following L^-boundedness of the vertical square function 





WfWp- 


B.2.2 Gradient estimates in Holder and As said above, we shall now work in the follow- 
Sobolev spaces ing setting, strictly weaker than the geometrical 

setting used in the first five sections of this work. 

Regularity assumptions 

(i) The metric measure space is doubling and the semigroup satisfies the 

Gaussian bound (lUEp . 

(ii) The gradient operator T satisfies (Gg^) and [DGq^) for some 2 < q < < oo, 

and the scale-invariant Poincare inequality. 

If go = 9 = 2, we require that the Davies-Gaffney estimates hold instead of 
{DGqp). As we shall see below, one can extend the machinery of paracontrolled calculus 
to that setting in Holder and Sobolev spaces. Recall the definition of the spaces IF 
and given in section im The parameter 6 is involved in the property {DGqfi). The 
following embedding is proved as ProDosition l2.5l bv using the fact proved in Lemma fB.il 
that the heat kernel is Holder continuous, with exponent 1 — 0, instead of its Lipschitz 
character. 


Proposition B.4. For a G (0,1), the space A*^ is continuously embedded into . If 
fj G (0,1 — 0), the two spaces are the same with equivalent norms. 

Sobolev spaces are naturally defined in terms of L as follows. 

Definition. Fix an exponent p G (l,oo), and s G M. A distribution / G 5', is said to 
belong to the inhomogeneous Sobolev space IT^’P if 


W^’P ■ = 


{1 + L)y ~||e-V|| + (1 + L)V 


< oo. 


► Proposition B.5. For a G (— oo, 1 — 0) and / G C®", we have 

sup (/ \ViTe~^^f\'^ dp] <t^||/||c-- 

x&M \ J/ 

The same conclusion holds with any of the operators , with an integer a > 1, in the 
role of 
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Proof — Consider 6 > 1, and write 


Vtr, 


-*^/= [ VtrQ^e-*^/ —+ 
Jo s 


For s < t, we have 

(-[ VtrgWe-^^/ 

V J B(x,\/i) 


dn 


(j^XU 

\S + tJ \JB(x,Vi) 


dfj, 


< 


)‘E^d/ 

QXtf 

iXo y2‘B{x,Vt) 

2 


dfi] , 


where are exponentially decreasing coefficients and where we used L^-L^ off- 
diagonal estimates of Fe””'^ (at the scale \/s + t ~ \/t) with the relation 


Q?i = 2''e-^V£t. 


So we have 


/ Ix/tiiQf'e-Vfrf,. <(f)‘E^'ll<3 


(b) 


~ \t) 

and we can integrate this inequality on the interval s € (0,t). For s > t, we use 
Property (DGg^e) to have 


/ VtTQXe-^^f'dS <(-] ^ (/ V^FQ^e-^V 

J B{x,y/t) / V'®/ \J B(x,y/s) 


d^ 




i-fl 

t\ 2 


Q{b+l)^-tLf 


L°°{B{x,^)) 


\ ^ 

<(^j s^WfWc-, 


where we have used qX = with L^-L^ (resp. L^-L°°) off-diagonal 

2 2 

estimates for (resp. {sL)qX/X)^ provided h is large enough. This inequality 

can be then integrated along s G (t, 1) as soon as 0 -|- cr < 1. 

We perform the same analysis for the term y/tVf, which gives 

1 

){b)-tLf\'l 


\ViTPXe-^^f\U^i] iFP^e-'^rd/i 

B{x,Vi) J \J B{x,l) , 


+ t 2 


LP^e-^^f 


L°°{B{x,l)) 


< t 2 


The conclusion follows from this inequality since t G (0,1) and a < 1 — 9. 


\> 
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Proposition B.6. For a S (0,1 — 0) and 0 < 25 < 1 — 0 — a, we have uniformly in a: G M 
and t > 0 


l-f \Vire-^^f\uS <M\{tL)^Qff 

\J B{x,y/i) J L 2 . 

where Ai is the Hardy-Littlewood maximal function. The same conclusion holds with any 
of the operators Pi°'\ with an integer a > 1 , in the role of and also by replacing 
B{x,y/i) by any bigger ball B{x,Ky/t) and the estimates are uniform with respect to 
K>1. 


Proof - We write (for a chosen large enough) 

roo 

\/fre-*^ = VtTL-'^e-^^L'If = Vi 

Jo 

For s < t, we then write 


ds 
s 2 




g(a-5)^_,Lg(5) 


and using off-diagonal estimates of the carre du champ of the semigroup, 

this yields 


(S \VirQVe-^^Lff\Ufi 

\J B{x,Vi) j 

For t < s we have by {DGqV 

1 

/ \VirQVe-^^L^f\Uf, 

J B{x,y/t) ; 


< 

rsj 


s\a 


M 


L^Qff] (x). 


< I - 

r\_; 


0-1 
s\ — 


V B(x,ys) 


-F 


. s-i 
•S \ 2 


< 

rsj 


5 \ 2 


^+<5 


M 


Q(a+l)^-tLL^j 

Ltgf/l (x), 


L°°{B{x,^)) 


where we used that with L^-L^ (resp. L^-L°°) off- 

diagonal estimates for Fgi“ (resp. provided a is large enough. Hence, 




VtTe-^^f 


dfJ-j < 

<M 


s\“ ds V /s\^+'5 ds 


L (VVfV (I) 


M 


L^QVf\{x) 


itLVQVf] {xl 


due to a large enough and 25 < 1 — 6 — a. We let the reader check the straight¬ 
forward modifications that are required to deal with a bigger ball B{x, Ky/i), and 
that the estimates are uniform with respect to iF > 1 . 

O 


Replacing the off-diagonal estimates by L^-L^ estimates, the same proof as 

above leads to the following result. 
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Proposition B.7. Assume the local Ahifors regularity. Suppose p G [l,oo) and a G 
(— 00 ,1 — 0 + -), and / G W^’^. Then, uniformly in x G M and t G (0,1], 


/ ViTe-^^f -dp] < 

K JBix.y/i) 




The same conclusion holds with any of the operators , with an integer o > 1, in the 
role of 

Proof — As previously, we write 


s 2 


for a a large enough integer. For s < t, we then have 


and for t < s we have by {DGq^) 


-/ y/tVQ^f^e dp\ < (y) t 2r||/||ry-. 

J B(x.\/i) / \t / 


/ ViTQi'^'^e-^^L^f'dpy <(j) ^ (/ iV^rg^e-^^Lf/l^d/x 

Q(^a+l)e-tLL§f 




L°°{B(x,^)) 




For the low frequency part, we have 


\.J B(x,\/i) 


1 . \ ^ , hzl, 

dp] < f 2 l|/|Hy«,p. 


Hence, 


[£ 


B{x,y/t) 
< 


VtTe dpj 

/ s \ O’ (f/s /*^/s 

i © (? 


0—1 u 7 

2 2p Cits . 1 -^ 


1 -- 
S 2 


Pi~ 




_ _i_ ^ 

< t 2 p + 2 


W“.P, 


due to a large enough and a < I — 9 + ^. 


t> 


--h| O' 
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B.3 Paraproduct and commutator estimates in This subsection is devoted to the state- 
Holder-Sobolev spaces ment and proofs of the main estimates 

about Paraproducts and commutators, 

in the current more general framework. 


B.3.1 Paraproduct estimates We state in this paragraph the basic continuity estimates 

satisfied by the maps defined by the low frequency part, the 
paraprodnct and the resonant terms - see Subsection 13.11 for the precise definition of 
these quantities and for detailed proofs. The low-frequency part is easily bounded. 

► Proposition B.8. Fix an integer b > 2. For any a,/3 S M and every 7 > 0 we have for 
every / S C" and 5 

(B.3) \\A_,{f,g)\\^^<\\f\\c^\\g\y. 

If the space (M, d, /r) is locally Alhfors regular, then for every a, /3 ,7 G M and p G [1, 00 ), 
we have for every / G W°‘’‘P and g 

(B.4) ||A_i(/,p)||^^,p < \\f\\w‘-,p\\g\\w^-p- 


The continuity properties of the paraproduct are given by the following statement. 


Proposition B.9. Fix an integer b >2. For any a G (—2,1) and / G we have 
• for every g G 


(B.5) 


(B.6) 


n(')(/) 


< 


lldllooll/lk 


for every g G with /? < 0 and a + {3 G (—2,1) 


n®(/) 


< 


Qa + P 


McM 


The proof is already given for Proposition 13.21 - and only relies on (jUED (which is 
also assumed here). We then state the analog in Sobolev spaces. 

Proposition B.IO. Assume local Alhfors regularity. Fix an integer b>2 and p G [l,oo). 
For any a G (—2,1) and / G we have 

• for every g G with ^ < /3 < 1 


(B.7) 


(B.8) 


p 

n®(/) 


< 




\\9\\wl^Af\\w°‘’P 


for every g G with /3 < ^ and a + P — ^ G (—2,1) 


n®(/) 


^a+/3-f,P ~ 


P 

WP’P 


lYa,p . 


Even if the proof is not very difficult, we give the details here in order to explain 
how to nse the L^-orthogonality property put forward in Lemma IB.31 


Proof - Recall that 
1 / 


(b) 


(q 




nf(/) = - / (tL)A 

lb Jo 

With s = a + {P — ^)_ > —2, Lemma [R3l yields 


.(b-i) 


((tL)pf / • Pf 5 


dt 

t 


l|n(')(/)|| 


W‘^’P ^ 


Q(b-l)j . 




+ t 




■f)' 
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(3 > v/p (and so s = a) then uniformly with respect to t > 0 we have due to 
the local Ahlfors regularity (which allows us to use a Sobolev embedding, see [T^ 
Lemma 10.5]) 


< lls-lloo < 


and so 

l|n(')(/)|| 




< 


gf-'V 




+ t 








IW^’P 


where we used again Lemma fB.31 If /3 < ^ (and so s = a + P — then 


d‘’slL< / l|ef'9lLv + llh‘’9 


ds 


y{b). 


< [ S^llgi'’ L2 g\\^— + \\g\\y^rp,p 
Jt S 


Hence, we conclude with Lemma lB.31 since 


l|n<‘>(/)|| 


< 


< 






+1 




2 2 
t 


I 


t> 


Proposition B.ll. Fix an integer b > 2. For any a, /3 £ (—oo, 1 — 0) with a + /3 > 0, for 
every / S C“ and 5 we have the continuity estimate 


n(')(/,5) 


< 


C°‘+P 




Proof — We only study the most difficult term in the resonant term n(^)(/, 5 ), which 
takes the form 

(B.9) A{f,g) := (VtPf "'V, VtPf "'^ 5 ) j- 

pj:^ satisfies L^-L°° off-diagonal estimates at order N {N can be chosen arbi¬ 
trarily large, since b is an integer) and so 


£>o 


\h\dp . 


l2^B{x,y/i) 

With h = y/iTP^^~^^ f ■ y/tTP^^~^^g and Holder’s inequality, we deduce that 


A{f,g){x 


e>0 \JB{x,2‘Vi) / \JB{x,2‘Vi) 


We then conclude as previously, with Proposition IB.51 


dt 


t> 
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We then give the analog estimate in Sobolev spaces. 

Proposition B.12. Assume the local Ahifors regularity. Fix an integer b > 2 and p G 
(l,oo). For any a, /3 G (—oo, 1- 
we have the continuity estimate 


(1, oo). For any a,l3 £ (—oo, 1 — 9) with a + /3 > ^, for every / G and g G W^’P, 






W<^’P llfflliv/3.p 


Proof — Again, we only study the most difficult term A{f,g) defined in (|B.9p . With 


s := a + /3 — ^ > 0, we have by Lemma fB.31 




< 




2 2 


Since s > 0, satisfies L}-L°° off-diagonal estimates at order | (see 

Proposition IB.2P and so 

\J 2^ Bix.'Jt') / 


With 


h = VtPPf 


and Holder’s inequality, we deduce that 
(B.IO) 


g) 


is bounded by 


E 2 

^>0 


—Is 


Jo \JB{x,2‘Vi) / \JB(x.2‘Vi) 


dt 


B{x,2^Vi) 


Then using Proposition IB. 71 with the Ahifors regularity, we have 

\JB(x,2‘Vt) J 

By Proposition IB. 61 we get 

(-f I ViTPt'^ /1' ) <M\itL)f Q® /] (x), 


lB(x,2‘Vi) 
for some (5 > 0. Hence 


g) 


< 

rs_/ 


< 


IW/5 


,p 

e>o 


M 




2 13—IL dt 
f P — 
t 


IWl^’P 


M 


Q^Ltf 


2 (iO “ 
t 


Using the Fefferman-Stein inequality on the maximal function together with the 
L^-boundedness of the square function - see Lemma fB.31 we deduce that 
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By a similar reasoning, we have 

rl 


\Mf, ff)||p ^ 


'0 


Vtpf Vipf "^^5) 


dt 
p t 




< 

rv-/ 


< 


< 




M 


Qf'Lif 


p dt 

t '2pt^ - 

p t 


M7“.p||s'llvi/^’P 

W<^’p\\9\\wf<-p^ 


'0 


^ s dt\ 

"T 


where we used that s > 0 and the L^-boundedness of the approximation operators. 
That concludes the proof of the estimate 

\\Mf, 9 )\\ws,P ^ WfWw^^’phWwP’P- 

Since the resonant part can be split into a finite number of terms similar to 
^{fj9)j we then deduce the Sobolev boundedness of the resonent part. 

O 


B.3.2 Commutator estimates We now focus on the adaptation of the commutator esti¬ 
mates given above in Proposition 13.41 

► Proposition B.13. Consider the a priori unbounded trilinear operator 

C{f,g,h) :=n('')(nf(/),/r) -^n^a/r), 

on Sg. Let a, (3 ,7 be Holder regularity exponents with a e (—1,1 — 0), /3 € (0,1 — 0) and 
7 e (— 00 , 1 ]. If 

0 < a + (3 + J and a - 1 - 7 < 0 

then, setting 5 := (a -|- /3) A (1 — 0) -|- 7 , we have 

(B.ll) \\C{f,9,h)\\cs<\\f\\c‘^\\9\\c^\\h\h, 

for every f ^ S and h S so the commutator defines a trilinear map from 

C^xC^ X a to C^. 

Proof — We refer to the proof of Proposition 13.41 for details and we keep the same 
notations. So it suffices to focus on a generic term of the form 

D{f,g,h) :=TZ{A{f,g),h) -gn{f,h) 

and prove the continuity estimate (jB.lip for it. As previously done, we split the 
proof of the commutator estimate (IB.Ill) for D in two steps, and introduce an 
intermediate quantity 

S{f,g,h) := (r{ViV^f,Vivfh)-Ptg) j, 

for which we shall prove that we have both 
(B.12) \\gn{f,h)-S{f,g,h)\\^, < \\f\\c<^ \\g\\cP ||h|b 

and 

(B.13) \\D{f,g,h)-S{f,g,h)\\^s < ||/lb \\h\\c^. 
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Step 1 — proof of ()B.12|1 . This part is very similar to Step 1 of Proposition I.S.41 
so we only point out the modifications. Using Gaussian pointwise estimates for the 
kernel of P/, we have for almost every x G M 




{ViVff,ViVfh)-{g{x)-Ptg)Yx) <^e 


£>0 


sup \g{x) - Ptg{y)\ 

yd{x,y)<2^Pt 


[i 


2^B(x,Vt) 


\ViT{V^f)fdfA l-f \ViT{V^h)\^dg] . 

J \J2‘B{x,Vt) J 

By using the C^-regularity of g as well as Proposition IB. 51 to estimate the 
averages, we get 


Vl(T{Viv'^f,ViVfh) ■ {g{x) - Ptg)Yx) < (4^t) 21212 j ||/||c“lbllc/3ll^llc^ 

Q: + /3 + 'y 

^ i^“ll/llc“||5'|lc/5|l^llo- 

Consequently, the continuity estimate (IB.12P in L°° comes from integrating with 
respect to time, taking into account the fact that a + /3 + 7 > 0. 

Then to estimate the regularity of gTZ{f, h) — S{f, g, h), one can exactly reproduce 
the same reasoning as for Proposition 13.41 bv using the Holder regularity of the heat 
kernel f Lemma IB. II) . which involves the condition a + /3 + 7 < 1 — 6 (since /3 < 1 — 6 
and a + 7 < 0 ). 

Step 2 — proof of ()B.13p . Given the collection := Qr^^)j.g(oi] operators, 
we need to prove that we have 


(B.14) 

for every r G (0,1], and where 
- S{f,g,h) 

(B.15) 


Qr(TZ{A{f,g),h) -S{f,g,h)^ 


^ 4 
^ r 2. 


= v^T (vt 1^' V^QI {Qlf • Ph) ^-Ptg- , ^tvfh 


dt 

t 


We are going to follow the same argument as for Proposition 13.41 and we only detail 
the modifications. So we set 

ds 

U 

t ■ 


Mf,g) := ViT V^Ql {Q^f • V^g) - - ngVff 


and using the L^-L°° off-diagonal estimates of P/, we then deduce that for almost 
every xq G M 


(B.16) 


'P{A{f,g),h) - S{f,g,h)^{xo) 

•-U’-'ii 


B(xo,2iPi) 


\At{f,g)\^dg 


'f 

\Je 


B(xoXVi) 


\ViTVfhfdg 


A, 


Using a suitable normalization of the operators QlQ^f ^ = / — Vif for some 
operator Pi, it yields for every x G M 

tif,g){x) < VlTV‘^Ql{Qlf{Vlg - Ptg{x))){x) ^ + \Ptg{x)\ViT{V^^V^f){x). 






















This quantity will then be integrated on := so we first aim to 

replace Ptg{x) by Ptg dfi. Observe that 

^ ds 


Mf,g)ix) < l^'Vtrv^Ql {Qlf{vlg-j^ Ptgdfi)"j (x) 


.17) 


+ 


ViTVffix) + \Ptgix)\VtT[V^Vif]ix). 


Ptg{x) - f Ptg d^i 

JBe 

As before, we use fi > Q and the regularity of g to have 

Me, 


Ptgix) - f Ptgdg. 

JBi 

and uniformly in y S M 


Plaiv) -i Ptgdfi < (max(s,t )2 +d{xo,y)^) \\g\\e. 

JB{xo,Vi) ^ ^ 

Moreover, it follows from the composition of off-diagonal estimates (correspond¬ 
ing to a analog of Lemma IA.5I - Parti, see also m Lemma 2.5]), that the 

operator y/tT(p^Ql^ satisfies off-diagonal estimates at the scale max(s,t) with 
extra factor ( max(s’t) ) ’ ^ ^8ix{s, t), we have with (jB.17j) 


an 


[L 


B{xo,2^Vi) 
1 


< 


\Mf^9)\^dy 
f mMsAy\ -ci’^tAk+e^ 


io ^„Uax(s,t)j ^ ^ [L 


ds 


IQifl^dy] Me 

B(xo,2^+^y/T) ) 'S 


+ 


[i 


< 


B(xo,2fVi) 

^ / min(s,t) 


\ViTP?Mdyj (2^Vi) ||5||c^+i^llffl|oo||/||c= 

13 ds 


S2t2 - \-t2 

s 


ll/llcHl5llc/^ 


/o \max(s,t)_ 

where we used Proposition IB.51 and the fact that a > —1 and a + j3 < 1 — 9 io 
estimate the integral over s. Observe that in the case where a-|-/?>l — 0, we get 


/ lAAf.gfdn] <2*’i^||/||c.||9llc». 

\JBixo,2^Vt} J 

Coming back to the identity ()B.16P with Proposition IB. 51 we have 

'■1 a + 0+'1 dt 


\\f\\cA\9\y\\h\\ 




P{A.{f,g),h) - S{f,g,h)J{xo) < P ^ 

A ll/l|c“||fl'llc/3|l^l|cT', 

since a -|- /3 -)- 7 > 0, uniformly for every xq G M. We then conclude to 

n{A{f,g),h)-s{f,g,h)) < WfMWaMWhM. 

Moreover, taking into account that we have qI^^PI = for t > r, we see 

that the estimate ()B.14p holds true (see thee proof of Proposition 13.4p . 


t> 
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We then aim to have a similar commutator estimate in Sobolev spaces. 

► Proposition B.14. Assume the local Ahifors regularity. Let a,/?, 7 be regularity exponents 

and p G (1, 00 ) with a G (—1,1 — 0), /3 G (^,1 — 0) and 7 G (— 00 ,1]. If 

2v V 

— < a + p + 7 and a + 7 < — 

P P 

then, setting (i:=(Q! + /3 — ^)A1 + 7 —^>0 and assume that 26 > (3 — We have 

(B-18) ||C'(/, 5, /i)||^,5_j, < ||/||vi/“.p Ibllw/^.p 

for every / G ,g G and /i G so the commutator defines a trilinear map 

from X W^’P x to W^’P. 

We follow the exact same proof as previously, so we keep the same notations and 
only focus on the modifications. 

Proof — Consider a generic term of the form 

D{f,g,h) :=n{A{f,g),h) -gn{f,h) 

and prove the continuity estimate (IB.lSp for it. Aiming that, we split into two 
terms by introducing the quantity 

S{f,g,h) := (r{VtV!f,VtVfh) • Ptg) j. 

for which we shall prove that we have both 
(B.19) \\gTZ{f,h) — S{f,g, /i)||^^_p < ||/||vv“.p ||5||vi//3,p ||/i||vyT.p 

and 

(B.20) \\D{f,g,h) — S{f,g,h)\\^sp < ||/||w“.p Hsllw/^.p II^IIv^t-.p- 

Step 1 — proof of (IB.19h . We first prove a weaker version of the continuity 
estimate (jB.19h . under the form of the inequality 

(B.21) WgPif, h) — S{f,g, h)\\^ < ||/||w“-p Ibllw/^’P II^IIwt-p- 

As previously, we have 

(B.22) (ff7^(/, h) - S{f, g, h)) {x) = Vj (T{VivU, Vivfh) • {g{x) - Ptg)) (x) j, 
for p-almost every x G M. Since g G , with /3 > ^ then g G W^’P C p and so 

\\Ptg-g\\oo < t 2 (^"p) \\g\\^,p,p. 

Hence 

/3 ^ 

(B.23) \Ptg{y) - g{x)\ < (yip d{x,y)) \\g\\wp,p, 

for every x,y G M. Coming back to equation (jB.22p and using Gaussian pointwise 
estimates for the kernel of Pi, we have for almost every x G M 

Vl{T{VtVlf,VtVlh) • {g{x) - Ptg)){x) 

is bounded above by 

f/ \VtT{Vlf)\^di^ (-f \VtT{vlh)\^dg] . 
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So using Propositions IB. 61 and IB. 51 we deduce that 


Vl Vivfh) ■ {g{x) - Ptg)){x) 


< M.[L2 f^[x) 11(711 vy/3.p II^11 


q: + ^+7) 

2 


P . 


Then the continuity estimate (IB.211) comes from integrating with respect to time, 
taking into account the fact that a + /3 + 7 > 


Let ns then estimate the regularity of F := gTZ{f, h) — S{f,g, h). It is known, see 
m Section 2.1.1], [21 Section 5.2] or m Proposition 9.7], that 


\\F\\ws,p<\\F\\p + \\Ss{F)\\p, 

where Ss is the Strichartz functional of index 6 £ (0,1): 

Ss{F){x) := I ( / \F{x)-F{y)\df,{y) 

\ -/O \JB{x,r) 

Fix r > 0 and two points x,y £ M with d{x, y) < r. Then as previously, we write 
F{x) -F{y) = [g{x)n{f,h) - S{f, g, h)){x) - [g{y)'R{f,h) - S{f,g,h){y)){y) 
=:U + V, 



with U defined by the formula 

{vl (riVlF^f, Vtv^h) ■ {g{x) - Ptg)) (x) - Vl VtV^h) ■ {g{y) - Fg)) (j/)} j, 

and V by 

r*! 

^ [vl (viViVlf, Virfh) ■ (gix) - Ptg)) (x) - V} [nViVlf, ViFfh) ■ {g{y) - Fg)) (j/)} 

By repeating previons arguments, we easily bonnd U as follows 


dt 

t 

dt 

t 


U < 


q: + /3+7 u 

t 2 p 


M{Ql^hff){x)+M{Qi^hff){y)\ Wgllw^AHw-r., 


for some e > 0 satisfying a + e < 1 — 9. For the second part, we nse the estimate 

|P| <a + b, 

with A equal to 


[vl (T{ViVff, Vivlh) ■ {g{x) - Fg)){x) - Vl (r{ViVlf, Vtvfh) ■ {g{x) - Fg)) ( 2 /)} 


and 


B := r \g{x)-g{y )\• \fI (v{ViFf f, Vivf h)) (y) 


dt 

t 


The last quantity is bounded (following the same estimates as previously since 
g G C^~p) by 

' f th[»+^-^)M[QfL^f]{y)- 

\J r2 


B < 


wP,p\W\w-^’P- 


For the quantity A, we combine the previous argument with the Holder regnlarity 
of the heat kernel, Lemma lB.il to get the npper bound 


A < 


rl / \ 1-e 

( / j. 


Vi 


||ff||H^.,p||/.|k..p. 
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The combination of all the previous estimates yields 
\F{x)-F{y)\ <\U\+A + B 


< 

rsj 


'j: 


tS/2 


M{QfL^^f){x)+M[Qf^L^^f]{y) 


(s) r # 


dt 


+ 


+ 


/ \ 1-0 
/ j_\ s 


Wt 


f){y) J 


dt 




l^^M[Qi^h^f]{y)j 


IIs'IIw/^.pII^IIvft'.p- 


This estimate holds uniformly for every y S B{x, r) and so can be averaged on this 
ball, 
that 


ball. We then conclude by Hardy’s inequality (with (5>0, (5<1 — 0 and /3 > ^) 


Ss{F) 


< 




2 2 
t 


W/s.pITIIVKt.p. 


Using Fefferman-Stein’s inequality and the L^-boundedness of the vertical square 
function (see Lemma fB.31) . we then deduce that 

ll-^llw'5.p ~ W^Wp + ll‘S'5(.P')||p ^ ||/||iV“.p|bllw/3.p||^l|iyT'.p, 
which concludes the proof of the continuity estimate (IB.191) . 

Step 2 — proof of (lB.20p . We need to prove that we have 


(B.24) 


* : = 


5 

L 2 


< 




{TZ{A{f,g),h)-Sif,g,h)) 
where TZ[A{f,g),h) —S{f,g,h) is equal to 

(B.25) V^T (Vt 1^' V?Ql {Qlf • Vlg) ^ - PtQ ■ , VtVth 


dt 

T' 

Using Lemma IB.31 with the L^-L°° off-diagonal estimates of V};, we deduce that 
quantity (*) is bounded above by a constant multiple of 


E 2 


-£5 


£>0 




B{xo,2‘Vi) 


\At{f, g)\‘^dy 


7 


dt 


B{xQ,2^Vt) 


\Virvfh\'^dy] — 


with 


At{f,g) :=Vir(^iy^Ql [Qlf-Vlg) ^ 


-Ptgvtf]- 


Following the reasoning in the previous proof of Proposition IB.131 by combining 
with Proposition IB.61 we can obtain that for almost every xq 


7 

\Je 


B{xo,2^A) 


Mf,9pdfi] < 2 ^^ p)tPd p)M[QP{tL)2 f]{xo)\\g\\^,p,^ 


+ t ^ M.[q['^ L'2 f]){xo)\\g\\wP’P 


Hence, since 25 > B — - we obtain 

’ r- p 






+ 
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We then conclude to (IB.24h by the Fefferman-Stein inequality and Lemma IB.31 


O 


B.3.3 Composition estimates The above continuity estimates are the main estimates used 

in the main part of this work to prove the paralinearisation 
and composition estimates for paraproduct. We state these results here in Holder and 
Sobolev spaces under our relaxed assumptions and leave the reader the task of checking 
that the proofs of section [3^ are easily adapted. 


Theorem B.15. Fix an integer b >2 and a nonlinearity F e Cl. 

(a) Let a S (0,1 — 0) be given. For every / S C", we have F(/) S C" and 

RfU) :=F(/)-n(f)^)(/)GC2“. 

More precisely, 


F(/)-ng)^^(/) 


< 


C2o 


|f|Ic3(i + 


If F G C^, then the remainder term Rpif) is Lipschitz with respect to /, in so far 
as we have 

||Ff(/) - RF{g)\\c2o. < l|F|b4 (1 + ll/lb + llffib)' 11/ - 5llc- 


(b) FixpG(l,oo). For every a G 1 — 0) and every / G we have F(/) G 


and 


Rpif) := F(/) - G 


2 


More precisely 

If F G C^, then the remainder term Rp{f) is Lipschitz with respect to /. 

Let us now examine the composition of two paraproducts. Note that for tt G C" and 
V € , with a G (0,1), /3 G (0, a], we have uv G 

About the composition of paraproducts, Theorem 13.61 in Holder spaces still holds 
since it only relies on the Gaussian estimate (lUEp : its Sobolev counterpart also holds. 
Theorem B.16. (a) Fix an integer 6 > 2, a G (0,1), (3 G (0,a] and consider u ^ 


and V G . Then for every / G we have 

nW(nW(/))-nW(/)GC“+^ 


with 


nW(nW(/))-nW(/) 


< 




lkllc“ 


(b) Fix an integer b > 2 and p G (l,oo). For a G (0,1) and (3 G (^jCrj. consider 
u G and v G W^’^. Then for every / G FF"’^, we have 

f«+/3-2^ ,P 


n(f)(nW(/))-nW(/)GW“ 





w°‘’P IpIIItIIvviS.p- 


with 
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B.3.4 Schauder estimates Proposition 13.81 gives an elementary proof in our setting of a 

Schauder-type estimate about the regularizing character of the 
convolution operation with the operators The same properties hold in our minimal 
setting since hey only rely on Gaussian property (lUEp and the semigroup structure, 
together with a Sobolev version which we state here without proof as it can be proved 
along the lines of proof of proposition 13.81 (Another approach can be also obtained 
by interpolating between the trivial case e = 1 and the limit case e = 0. The latter 
case e = 0, corresponds exactly to the so-called maximal regularity which has been 
the topic of a huge literature, see for example [33] where the Gaussian upper estimates 
(lUEh are used.) 


► Proposition B.17. (a) Consider /3 e M and e € (0,1). For every T > 0 and v € 

CtC^, the function V{t) := P^’^gV{s) ds belongs to with 

||E(t)||^3+2_2, < sup ||?^(s)||c /9 
se[o,r] 


and 


So 


||^|IcW+2-2£)/2^<^ ~ II^IICtC/3- 


<rn|/|b^C«- 2 + 2 e. 


(b) Consider /3 e M, p e (1, oo) and e G (0,1). For every T > 0 and v G the 

function V{t) := P^^gV{s) ds belongs to with 

||^W|lw^+2-2..p < r" sup ||?^(s)||^y,,,p 
se[o,r] 


and 


^11,^1 (/9+2-2e)^^ 




B.4 Resolution of PAM equation in such a Building on the estimates proved in this 
2-dimensional setting Appendix, it is elementary to introduce and 

study paracontrolled distributions in Holder 
and Sobolev spaces along the lines of Subsections 14.11 and 14.21 in the present extended 
setting. Its application to the parabolic Anderson model equation (PAM) is also almost 
straightforward as we only need to check that the renormalization procedure explained 
in details in subsection 15.31 under the (Lip) assumption can be run here as well. This 
is indeed the case if the exponent qo in the gradient assumption (G^g) is large enough, 
as this assumption yields some -Gaussian” estimates for every qi < go- 

Let us compute, as an example, an integral of type 


Is,t-= / T^pt{x,y)r^ps{x,z)dn{x) 
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where pt is the heat kernel of e and s < t. By (IB.lh with the local Ahlfors regularity, 
there exists a positive constant c such that we have 




/ 


d(x,y)^ d(x,z)^ 

p{dx] 


1 -^ 

n 


<t 291 / Qt{x,y)Qs{x,z) p{dx) 


1 — 


91 


< 


t 2915 n) ^ 


where we used Lemma fA.51 So with respect to Subsection l5..S[ where (Lip) was assumed 
and where Is,t would be estimated by gs+t{y,z), we now have the estimate 
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Gt+s{y,z), 


involving an extra factor • Since all the conditions on the exponents were 

open conditions in Subsection 15.31 we may allow a small loss if it is small enough. As a 
consequence, we deduce that if qi can be chosen large enough then we may adapt and 
repeat the renormalization procedure of the white noise in Holder and Sobolev spaces. 
The latter condition on qi is equivalent to taking qo big enough. 


We summarize this result under the following form, which gives an analogue of 
theorem O 


► Theorem B.18. Assume the local Ahlfors regularity of dimension 2, as well as {P 2 ), (Gg^) 
and (DGqp) for <70 lai'ge enough and 6 small enough. Fix p > 2 a large enough exponent. 
Let ^ stand for a time-independent weighted noise in space, and set := and 
= ds. 

(a) The pair converges in probability in some space (in the Holder scaling 

(C^)s or Sobolev scaling to some extended noise (C^Al), with C = C and 

n(A, C) well-defined in the above sense. 

(b) Furthermore, if stands for the solution of the renormalized equation 

(B.26) dtu^ + Lu^ = — P F'[u^) F{u^), u^{0)=uo 


where P{-) 


E 




is a deterministic real-valued function on M, 


then converges in probability to the solution u of {gPAM) associated with 
(C,A1), in some space whose definition depends on whether or not F is linear. 


References 

[ 1 ] P. Auscher, Regularity theorems and heat kernel for elliptic operators, J. London Math. Soc., 
54 (1996), no. 2, 284-296. 

[2] P. Auscher, T. Coulhon, X. T. Duong, and S. Hofmann, Riesz transform on manifolds and heat 
kernel regularity, Ann. Sci. Ecole Norm. Sup., 37 (2004), no. 4, 911-957. 

[3] D. Albrecht, X. T. Duong and A. McIntosh, Operator theory and harmonic analysis, in Instruc¬ 
tional Workshop on Analysis and Geometry, Part III (Canberra, 1995), Proc. Centre Math. 
Appl. Austral. Nat. Univ. 34 (1996), 77-136. 

[4] P. Auscher, A. McIntosh and P. Tchamitchian, Noyau de la chaleur d’opfateurs elliptiques 
complexes. Math. Research Lett. 1 (1994), 37-45. 













81 


[5] N. Badr, F. Bernicot, and E. Russ, Algebra properties for Sobolev spaces-applications to semi- 
linear PDFs on manifolds, J. Anal. Math., 118, no.2 (2012), 509-544. 

[6] I. Bailleul, A flow-based approach to rough differential equations, https://perso.univ- 
rennesl.fr/ismael.bailleul / files/M2Course.pdf 

[7] I. Bailleul, F. Bernicot and D. Frey, Higher order paracontrolled calculus, 3d-PAM and multi¬ 
plicative Burgers equations, arXiv:1506.08773, 2015. 

[8] H. Bahouri, R. Danchin and J.Y. Chemin, Fourier analysis and nonlinear partial differential 
equations, Grundlehren der Mathematischen Wissenschaften 343, (2011). 

[9] F. Baudoin and N. Garofalo, Curvature-dimension inequalities and Ricci lower bounds for sub- 
Riemannian manifolds with transverse symmetries, J. Eur. Math. Soc. (2015). 

[10] F. Bernicot, A T(l)-Theorem in relation to a semigroup of operators and applications to new 
paraproducts, Trans. Amer. Math. Soc. 364 (2012), 257-294. 

[11] F. Bernicot, T. Coulhon and D. Frey, Gaussian heat kernel bounds through elliptic Moser 
iteration, J. Math. Pures Appl. (2016), arXiv;1407.3906. 

[12] F. Bernicot, T. Coulhon and D. Frey, Sobolev algebra through Heat semigroup, submitted, 
arXiv:1505.01442. 

[13] F. Bernicot and D. Frey, Riesz transforms through reverse Holder and Poincare inequalities, 
submitted, arXiv: 1503.02508. 

[14] F. Bernicot and Y. Sire, Propagation of low regularity for solutions of nonlinear PDEs on a 
Riemannian manifold with a sub-Laplacian structure, Ann. I. H. Poincare - AN 30 (2013), 
935-958. 

[15] J.M. Bony, Calcul symbolique et propagation des singularies pour les equations aux derivees 
partielles non lineaires, Ann. Sci. Eco. Norm. Sup. 114 (1981), 209-246. 

[16] S. Boutayeb, T. Coulhon and A. Sikora, A new approach to pointwise heat kernel upper bounds 
on doubling metric measure spaces, Adv. in Math. 270 (2015), 302-374. 

[17] H-Q. Bui, X. T. Duong and L. Yan, Calderon reproducing formulas and new Besov spaces 
associated with operators, Adv. in Math. 229 (2012), 2449-2502. 

[18] R. Catellier and K. Chouk, Paracontrolled distributions and the 3-dimensional stochastic quan¬ 
tization equation. To appear in Ann. Probab., 2016. 

[19] T. Coulhon, E. Russ, and V. Tardivel-Nachef, Sobolev algebras on Lie groups and Riemannian 
manifolds, Amer. J. of Math., 123 (2001), 283-342. 

[20] T. Coulhon and A. Sikora, Gaussian heat kernel bounds via Phragmen-Lindelof theorem, Proc. 
London Math. Soc. 96 (2008), 507-544. 

[21] M. Cowling, I. Doust, A. McIntosh and A. Yagi, Banach space operators with a bounded 77°° 
functional calculus, J. Austral. Math. Soc. Ser. A 60 (1996), 51-89. 

[22] E. De Giorgi, Sulla differenziabilita de analiticita delle estremali degli integral! multipli regolari, 
Mem. Accad. Sci. Torino Cl. Sci. Pis. Mat. Nat., 3 (1957), no. 3, 25-43. 

[23] X.-T. Duong and D. W. Robinson, Semigroup kernels, Poisson bounds, and holomorphic func¬ 
tional calculus, J. Funct. Anal. 142(1) (1996), 89-128. 

[24] P. Friz and M. Hairer, A course on rough paths (with an introduction to regularity stuctures). 
Springer, 2014. 

[25] P. Friz and N. Victoir, Multidimensional stochastic processes as rough paths; Theory and Ap¬ 
plications. Cambridge University Press, 2010. 

[26] M. Fukushima, Y. Oshima and M. Takeda, Dirichlet forms and symmetric Markov processes. 
De Gruyter Studies in Mathematics 19, Walter de Gruyter, Berlin, 1994. 

[27] I. Gallagher and Y. Sire, Besov algebras on Lie groups of polynomial growth, Studia Math. 212 
(2012), no. 2, 119-139. 

[28] A. Grigor’yan, Gaussian upper bounds for the heat kernel on arbitrary manifolds, J. Diff. Geom. 
45 (1997), 33-52. 

[29] A. Grigor’yan, Analytic and geometric background of recurrence and non-explosion of the Brow¬ 
nian motion on Riemannian manifolds. Bull. Amer. Math. Soc. 36 (1999), 135-249. 

[30] A. Grigor’yan and L. Liu, Heat kernel and Lipschitz-Besov spaces. Forum Math. (2015). 

[31] M. Gubinelli, Controlling rough paths, J. Funct. Anal. 2016, no. 1, (2004), 86-140. 

[32] M. Gubinelli, P. Imkeller and N. Perkowski, Paracontrolled distributions and singular PDEs, 
Preprint 2012, arXiv:1210.2684. 

[33] Y. Guivar’ch, Croissance polynomiale et periodes des functions harmoniques. Bull. Soc. Math. 
France 101 (1973), 333-379. 

[34] P. Gyrya and L. Saloff-Coste, Neumann and Dirichlet heat kernels in inner uniform domains, 
Asterisque 33 (2011), Soc. Math. France. 


82 


[35] M. Hairer, A theory of regularity structures, Invent. Math., 198 (2014), no. 2, 269-504. 

[36] M. Hairer, Rough Stochastic PDEs, Commun. Pure Appl. Math., 64 (2011), no. 11, 1547-1585. 

[37] M. Hairer, Solving the KPZ equation , Annals of Maths, 178 (2013), no. 2, 559-664. 

[38] M. Hairer and C. Labbe, A simple construction of the continuum parabolic Anderson model on 
R^, Elect. J. Probah., 2016. 

[39] M. Hairer and C. Labbe, Multiplicative stochastic heat equations on the whole space. Preprint, 
arXiviarXiv: 1504.07162. 

[40] M. Hairer and H. Weber, Rough Burgers-like equations with multiplicative noise, Probab. Theory 
Pel. Fields, 155 (2013), no. 1, 71-126. 

[41] W. Hebisch and L. Saloff-Coste, Gaussian estimates for Markov chains and random walks on 
groups, Ann. Probab. 21 (1993), no. 2, 673-709. 

[42] W. Hebisch and L. Saloff-Coste, On the relation between elliptic and parabolic Harnack in¬ 
equalities, Ann. Inst. Fourier 51 (2001), no. 5, 1437-1481. 

[43] M. Hieber and J. Priiss, Heat kernels and maximal estimates for parabolic evolution 

equations. Comm, in Partial Differential Equation 22 (1997), 1647-1669. 

[44] S. Ishiwata, A Berry-Esseen type theorem on nilpotent covering graphs, Canad. J. Math. 56 
(2004), no. 5, 963-982. 

[45] S. Ishiwata, Gradient estimate of the heat kernel on modified graphs. Pot. Analysis 27 (2007), 
no. 4, 335-351. 

[46] P. C. Kunstmann and L. Weis, Maximal Lp regularity for parabolic problems, Fourier multiplier 
theorems and -functional calculus, Lect. Notes in Math. 1855. Springer-Verlag (2004). 

[47] H-Q Li, La transformation de Riesz sur les varietes coniques, J. Funct. Anal. 168 (1999), 
145-238. 

[48] P. Li and S. T. Yau, On the parabolic kernel of the Schrodinger operator, Acta Math. 156 
(1986), 153-201. 

[49] L. Liu, D. Yang and W. Yuan, Besov-type and Triebel-Lizorkin-type spaces associated with 
Heat kernels, Preprint 2013, arXiv:1309.1366. 

[50] T. Lyons, Differential equations driven by rough signals, Rev.Mat. Iberoamericana 14(2), 1998. 

[51] M. Caruana, T. Levy and T. Lyons, Lect. Notes Math. 1908, 2007. 

[52] Y. Meyer, Wavelets and operators, Cambridge Studies in advanced mathematics (1992). 

[53] A. Nagel, E. M. Stein and S. Wainger, Balls and metrics defined by vector fields I: Basic 
properties, Acta Math. 155 (1985), 103-147. 

[54] B. Qian, Hamilton type gradient estimate for the sub-elliptic operators. Potential Anal. (2015). 

[55] L. Saloff-Coste, Analyse sur les groupes de Lie a croissance polynomiale, Ark. Mat. 28 (1990), 
315-331. 

[56] L. Saloff-Coste, A note on Poincare, Sobolev, and Harnack inequalities, Int. Math. Res. Not. 
(1992), no. 2, 27-38. 

[57] E. M. Stein, Topics in harmonic analysis related to the Littlewood-Paley theory. Princeton 
University Press, 1970. 

[58] K.T. Sturm, Analysis on local Dirichlet spaces 1. Recurrence, conservativeness and L^-Liouville 
property, J. Reine Angew. Math. 456 (1994), 173-196. 

[59] K.T. Sturm, Analysis on local Dirichlet spaces H. Upper Gaussian estimates for the fundamental 
solutions of parabolic equations, Osaka J. Math. 32 (1995), no. 2, 275-312. 

[60] N. Th. Varopoulos, Analysis on Lie groups, J. Funct. Anal. 76 (1988), 346-410. 

[61] N. Th. Varopoulos, Small time Gaussian estimates of heat diffusion kernels. Part 1. The semi¬ 
group technique, Bull. Sci. Math. Ser. 2 113 (1989), 253-277. 

[62] N. Th. Varopoulos, L. Saloff-Coste and T. Coulhon, Analysis and geometry on groups. Cam¬ 
bridge University Press, 1992. 

[63] F.-Y. Wang and L. Yan, Gradient estimate on convex domains and applications, Proc. Amer. 
Math. Soc. 141 (2013), 1067-1081. 

[64] R. Zhu and X. Zhu, Three dimensional Navier-Stockes equation driven by space-time white 
noise, arXiv:1406.0047, (2014). 

[65] R. Zhu and X. Zhu, Approximating three-dimensional Navier-Stokes equations driven by space- 
time white noise, arXiv:1409.4864 , (2014). 


83 


Institut de Recherche Mathematiques de Rennes, 263 Avenue du General Leclerc, 
35042 Rennes, France 

E-mail address’, ismael.bailleul@uiiiv-rennesl.fr 

CNRS - Universite de Nantes, Laboratoire de Mathematiques Jean Leray, 2, Rue de la 
Houssiniere 44322 Nantes Cedex 03, France 

E-mail address’, frederic.bernicot@univ-nantes.fr 

CNRS - Universite Paris-Sud, Laboratoire de Mathematiques, UMR 8628, 91405 Orsay, 
France 

E-mail address: dorothee.frey@univ-nantes.fr 


